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ABSTRACT 

A literature survey on thermal stresses during welding and buckling 
after welding is described. From the analyses found, a computer program was 
developed to calculate thermal stresses and resulting residual stresses due to 
a moving heat source. This program is included in the appendix. 


Studies at Marshall Space Flight Center on three-dimensional movement 
during welding of flat plate specimens are summarized and interpreted on the 
basis of the Battelle literature survey and computer analysis. It is concluded 
that local metal movement during welding appears to be caused by bending 
moment and that general metal movement is apparently caused by buckling. 
Future programs to improve the present analysis are recommended. 
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FOREWORD 


This is the final report on a program to analyze weld distortion in 
Saturn V structures under Contract No. DA-0 1-02 1-AMC- 14693 (Z) for the 
Manufacturing Engineering Laboratory of the George C. Marshall Space Flight 
Center, Huntsville, Alabama. 

Information in this report comes from documents and data supplied by 
the Marshall Space Flight Center of the National Aeronautics and Space Admin- 
istration and presently available literature. 
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Section I. INTRODUCTION 


A major problem in the fabrication of space vehicle structures is 
distortion caused by welding. These structures include a number of structural 
components that must be fabricated to close dimensional tolerances. Unaccept- 
able distortion that occurs during welding of joints in these structures is often 
impossible to correct without removing the weld joint. Distortion can also cause 
mismatch of weld joints , which results in reduction of the joint strength. 

Empirical approaches to the problem of distortion have neither success- 
fully prevented distortion nor given an understanding of thermal and mechanical 
causes of distortion. Yet, before this program no organized analytically 
oriented research had been conducted on distortion control of space vehicle 
structures. 

Several aerospace companies have encountered distortion problems during 
fabrication of Saturn V components. Although production practices have been 
developed to overcome these problems, these again are empirical solutions. 
There has been no opportunity to develop theoretical knowledge of these 
problems. The objective of the Battelle study was to theoretically analyze 
certain distortion phenomena encountered during the fabrication of Saturn V 
structures and to correlate this analysis with empirical data developed during 
their production. The emphasis of this study was placed on the mathematical 
analysis of thermal stresses and on metal movement in the thickness direction 
during fusion welding. The study involved the following phases: 

1) Identification of weld distortion problems 

2) Literature survey and analysis of present data on the basis of 
presently available formulas 

3) Formulation of experimental programs to improve the theoretical 
analysis 

4) Recommendation of future programs on distortion control. 

The objective of the first phase was to identify weld distortion problems 
in the Saturn V structures. During discussions with National Aeronautics and 
Space Administration (NASA) monitors on this research program, it was 
decided that the Battelle study would be focused on the analysis of metal 
movement around the are during welding. 
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In the previous fabrication of Saturn V fuel tanks, mismatch of joints 
was experienced during the welding of girth joints, especially the joint between 
the cylindrical parallel section and the transition ring. Before the initiation of 
the Battelle study, an experiment was conducted at Marshall Space Flight 
Center (MSFC) to measure metal movement during girth welding of a 
cylinder 133 inches in diameter and 12 inches high. In this experiment it was 
found that there is metal movement in the thickness direction near the arc. A 
buckling distortion upon completion of the weld was also observed. Therefore 
it was decided to direct the Battelle research toward a study of these two 
problems, especially the metal movement during welding. 

A literature survey was made to obtain useful information on thermal 
stresses during welding and buckling after welding. Section III presents the 
results of this survey as a state-of-the-art analysis of weld distortion. Follow- 
ing the survey, the existing distortion data on the Saturn V structure were 
analyzed on the basis of presently available analyses. Using these analyses, 
a computer program for analyses of thermal stress was developed. Section IV 
describes results of the computer analysis. The FORTRAN computer program 
is given in the appendix. 


On the basis of the above survey, Battelle recommended that MSFC 
conduct studies on three-dimensional movement during welding of flat-plate 
specimens. Flat plates were recommended because they yield the most funda- 
mental information on thermal stress during welding and buckling after welding. 
Section V of this report presents summaries of the experimental results. 

Section VI discusses how the experimental results can be interpreted on the 
basis of the literature survey and the analysis conducted at Battelle. 

The entire program is summarized and pertinent conclusions are 
presented in Section VII. On the basis of the information obtained in this 
research, recommendations are made for future research programs. This dis- 
cussion is given in Section VIII, 
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Section li. GENERAL DISCUSSION OF FABRICATION OF SATURN V 


The Saturn V is a giant three-stage rocket designed to send the Apollo 
spacecraft toward the moon. It stands 362 feet high, weighs 6. 1 million 
pounds , and has a maximum diameter of 33 feet. Its fabrication includes 
thousands of feet of welding and introduces many welding problems. Several 
articles have been written describing the fabrication and welding of the Saturn 
V structure. 1-5 The following is a very limited description of steps in the 
fabrication of the first stage (S-1C) that are most directly related to the 
Battelle study. 

The S- 1C is a semi-monocoque structure measuring 138 feet in length 
by 33 feet in diameter. The structure has various components including welded 
fuel and oxidizer tanks. These tanks are fabricated from two aluminum alloys, 
2219 for the S-1C and 2014 for the S— II and S-IV stages. 


The Battelle study was related primarily to these fuel and oxidizer tanks 
which are of all-welded construction. Either gas metal-arc (GMA) or gas 
tungsten-arc (GTA) welding is used on these assemblies except for the inside 
baffles . Both tanks are constructed by attaching skin cylinders to upper and 
lower head assemblies by means of a Y-shaped transition ring. The unit serves 


as a transition piece between bulkheads 
ture and forward skirt. This design is 



B-B 


A-A COMPLETED WELD 

B-B STRAP ALIGNMENT CLAMP 


FIGURE 1. WELDING OF THE 
TRANSITION RING 


and skins and between the thrust struc- 
illustrated in Figure 1. 

Welding between the transition 
ring and cylindrical sections is probably 
the most critical on the S-1C stage and 
thus the most troublesome from a 
distortion viewpoint. Three vertical 
joints, 47 inches long, are automatically 
GMA welded to join three forged and 
rolled billets. An excess of 0. 140 inch 
of metal is allowed for weld shrinkage 
at each joint. During final assembly of 
the tanks, the transition ring is joined 
to the cylindrical parallel section by a 
GMA or GTA weld made from both sides 
in the horizontal position. 
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Section III. PRESENT STATE OF THE ART 
OF THE ANALYSIS OF WELD DISTORTION 


This section describes the state of the art of the analysis of weld 
distortion. Emphasis is placed on distortions that occur during butt welding of 
girth joints of a cylindrical tank. 


1. Past Reviews on Weld Distortion and Residual Stresses 

Residual stresses and distortion during fabrication of welded 
structures have been studied by many investigators, Although most studies 
have been experimental, some analytical studies have been made. Several 
reviews and books have been written on residual stresses and distortion in 
weldments and welded structures . 

Spraragen and associates prepared a series of comprehensive reviews 
on residual stress and distortion in welds in 1937 through 1950. 6-8 In 1954, 
Osgood 9 edited a book on residual stresses in metals and metal construction. 
Masubuchi wrote a report on measurement of residual stresses in metals and 
metal structures in 1965. 10 His report emphasized residual stresses in 
structures in aluminum alloys. 

Four books have been published in Japan on residual stresses and dis- 
tortion in weldments by Naka, Kihara, Masubuchi, Watanabe, and Satoh. 11-14 
In Russia several books have been written including one by Okerblom. 15 The 
Swede, Gunnert, has written a book that covers primarily the residual-stress 
measuring technique he developed. 16 A group of international authors, Hall, 
Kihara, Soete, and Wells, 17 published a book in 1967 that summarizes recent 
work on the effects of residual stresses on brittle fracture of welded steel 
plates. 


Among these past reviews, those which were written by Spraragen 
et al, 6-8 Kihara and Masubuchi, 12 > 13 Watanabe and Satoh, 14 and Okerblom 15 con- 
tain considerable information on the analysis of weld distortion. However, most 
of the above reviews cover only welding of mild steel with covered electrodes. 
Very little is mentioned about other materials, such as high-strength ferrous 
and nonferrous alloys, or inert-gas metal or tungsten-arc welding processes. 

Since 1 July 1967 Battelle has been conducting a 3-year program, 
supported by a group of industrial companies , to develop analytical methods for 
prediction of distortion in welded structures. 18 As a part of the program, 
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Battelle has recently compiled a very extensive and up-to-date bibliography on 
the subject of weld distortion. Battelle' s review covers many materials and 
many processes and includes literature from eighteen different countries. This is 
probably the most extensive collection of references available on the subject of 
weld distortion. Direct access to this review is restricted to supporting 
companies. 


a. General Discussion of Residual Stresses 
and Distortion in Welded Structures 


Residual stresses are those stresses that would exist in a body 
if all external loads were removed. Residual stresses in metal structures 
occur for many reasons during various manufacturing stages, including rolling, 
casting, machining, flame cutting, and welding. 

Because a weldment is locally heated by the welding heat source, the 
temperature distribution in the weldment is not uniform and changes as welding 
progresses. During the welding cycle, complex strains occur in the weld 
metal and base metal regions near the weld, both during heating and cooling. 
The strains produced during heating are accompanied by plastic upsetting. 

The stresses resulting from these strains combine and react to produce 
internal forces that cause bending, buckling, and rotation. It is these displace- 
ments that are called distortion. 


b. Typical Distributions of Residual Stresses in Weldments 

The residual stresses in a welded joint are caused by the 
contraction of the weld metal and the plastic deformation produced in the base 
metal region near the weld. Residual stresses in a welded joint are classified 
as "residual welding stress, " which occurs in a joint free from any external 
constraint, and "reaction stress'' (or "locked- in stress") , which is induced by 
an external constraint. Most of the information on residual stresses has been 
obtained on arc-welded joints, especially on butt joints in mild steel. Rather 
limited information is available on residual stresses and distortions in weld- 
ments in aluminum alloys. 19-22 

A typical distribution of residual stresses in a butt weld is shown in 
Figure 2. The stresses of concern are those parallel to the weld direction, 
designated o^, and those transverse to the weld, designated a . 


5 



Y 



Y 

fa) BUTT WELD 


Y 



Y 

(b) DISTRIBUTION OF a ALONG YY 



(c) DISTRIBUTION OF a ALONG XX 

y 


FIGURE 2. TYPICAL DISTRIBUTION OF RESIDUAL STRESSES 
IN BUTT WELD 
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The distribution of the residual stress along a line transverse to the 

weld (YY) is shown in Figure 2 (b) . Tensile stresses of high magnitude are 
produced in the region of the weld; these taper off rapidly and become com- 
pressive at a distance of several times the width of the weld. The weld metal 
and heat-affected zone try to shrink in the direction of the weld, and the 
adjacent plate material prevents this shrinkage. The distribution of a 

y 

residual stress along the length of the weld XX is shown by Curve 1 in Figure 
2(c) . Tensile stresses of relatively low magnitude are produced in the middle 
of the joint, and compressive stresses are observed at the end of the joint. 

If the contraction of the joint is restrained by an external constraint, 
the distribution of is as shown by Curve 2 in Figure 2(c) . Tensile stresses, 

approximately uniform along the weld, are added as the reaction stress. An 
external constraint, however, has little influence on the distribution of a 

x 

residual stresses. 


The magnitude and distribution of residual stresses in a weld are 
determined by ( 1) expansion and contraction characteristics of the base metal 
and the weld metal during the welding thermal cycle and (2) temperature versus 
yield-strength relationships of the base metal and the weld metal. Much 
research on carbon- steel weldments has shown that the maximum residual 
weld stress is as high as the yield stress of the weld metal. However, in a 
recent investigation conducted at Battelle, 23 the maximum stresses in high- 
strength- steel weldments were not as high as the yield strength of the weld 
metals. In weldments in heat-treated SAE 4340 steel, for example, the maxi- 
mum residual stresses were around 50,000 to 80, 000 pounds per square inch, 
considerably less than the yield strength of the weld metal (around 150, 000 psi) 
and the base metal (224, 000 psi) . 

Hill 20 investigated residual stresses in butt joints in Alloy 5456-H321 
plates welded by the inert-gas-shielded arc-welding process with Alloy 5556 
consumable electrodes. Figure 3 shows a typical distribution of longitudinal 
residual stresses in a 0. 5- by 36- by 48-inch panel fabricated by welding two 
0.5- by 18- by 48-inch plates. Variation in yield strength of the material with 
distance from the center line of weld is also shown. The residual tensile 
stresses in and adjacent to the weld approach the yield strength. These tensile 
stresses are confined to the region in which the heat of welding has lowered the 
yield strength of the material. 
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FIGURE 3. DISTRIBUTION OF YIELD STRENGTH AND LONGITUDINAL 
RESIDUAL STRESSES IN A WELDED 5456-H32 1 PLATE 


c. Fundamental Types of Distortion 


The distortion found in fabricated structures is caused by 
three fundamental dimensional changes that occur during welding, (1) transverse 
shrinkage that occurs perpendicular to the weld line, (2) longitudinal shrinkage 
that occurs parallel to the weld line, and (3) an angular change that consists of 
rotation around the weld line. All three of these dimensional changes are shown 
in the sketches in Figure 4. 

Figure 4(a) shows transverse shrinkage in a simple butt weld. The 
distribution of the longitudinal residual stress cr^ is shown in Figure 4(b) . This 

stress causes the longitudinal shrinkage shown in the figure. High tensile 
stresses exist in regions near the weld, and compressive stresses exist in 
regions away from the weld. This indicates that the longitudinal shrinkage of 
the weld metal and base metal regions adjacent to the weld are restrained by the 
surrounding base metal regions . 
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RESIDUAL STRESS 



(a) TRANSVERSE SHRINKAGE (b) LONGITUDINAL SHRINKAGE 
IN A BUTT WELD IN A BUTT WELD (DISTRIBUTION 

OF LONGITUDINAL RESIDUAL 
STRESS <T IS ALSO SHOWN) 



(c) ANGULAR CHANGE IN 
A BUTT WELD 


(d) ANGULAR CHANGE IN A FILLET 
WELD 



FIGURE 4. FUNDAMENTAL DIMENSIONAL CHANGES THAT OCCUR 
IN WELDMENTS 


Figure 4(c) shows the angular change that occurs in a butt weld. Non- 
uniformity of transverse shrinkage in the thickness direction is the cause of 
this rotation. Figure 4(d) shows the angular change that occurs in a fillet weld. 
Here, the angular change is caused by the unbalance of shrinkage on opposite 
sides of the flange member. 

Distortions which occur in practical weldments are more complex than 
those shown in Figure 4. Figure 5 shows distortion induced by longitudinal 
shrinkage in various joints. Figure 5(a) shows distortion of a beam fabricated 
by welding. When the weld line is not located along the neutral axis of the 
joint, longitudinal shrinkage of the weld causes bending distortion. Figure 5(b) 
shows the distortion of a single-vee butt weld. The amount of weld metal, 
and consequently the amount of shrinkage, is not uniform in the thickness 
direction. The weldment deforms into a concave shape (looking from the top) . 
When a thin plate is welded, residual compressive stresses, which occur in 
areas away from the weld, cause the buckling shown in Figure 5(c) . In the 
cases shown in Figures 5(a) , (b) , and (c) , distortions caused by angular 
changes are superimposed on those caused by longitudinal shrinkage. 
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(a) LONGITUDINAL DISTORTION OF A BUILT-UP BEAM 



(b) LONGITUDINAL BENDING DISTORTION OF A SINGLE-VEE . 
BUTT WELD 



(c) BUCKLING DISTORTION 


FIGURE 5. DISTORTIONS INDUCED BY LONGITUDINAL SHRINKAGE 


Figures 4 and 5 show only a few examples of weld distortion. Shrinkage 
and distortion that occur during the fabrication of actual structures are far 
more complex than those shown. A recent paper by Masubuchi 18 covers various 
types of weld distortion. 


2. Mathematical Analysis of Thermal Stresses 
During Welding and of Residual Stresses 

a. Changes of Temperature and Stress During Welding 


Figure 6 shows schematically how residual stresses are formed 
in a weld. Figure 6(a) shows a bead-on-plate weld in which a weld bead is 
being laid at a speed v. The coordinate axis is 0-xy, the origin 0 is on the 
surface underneath the welding arc, and the x direction lies in the direction of 
welding. 

Figure 6 shows temperature distribution along several cross sections. 
Along Section A-A, which is ahead of the welding arc, the temperature change 
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(4) SECTION D-D 


(b) TEMPERATURE CHANGE (c) STRESS 


FIGURE 6. SCHEMATIC REPRESENTATION OF CHANGES OF TEMPERATURE 
AND STRESSES DURING WELDING 

due to welding, AT, is almost zero. Along Section B-B, which crosses the 
welding arc, the temperature distribution is very steep. Along Section C-C, 
which is some distance behind the welding arc, the distribution of temperature 
change is shown in Figure 6(b) . Along Section D-D, which is very far from the 
welding arc, the temperature change due to welding again diminishes. 
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Figure 6(c) shows the distribution of stresses along these sections in 

the x direction, a . Stress in the y direction a and shearing stress t also 
x y xy 

exist in a two-dimensional stress field [ Figure 6 (a) ] . * 

Along Section A-A [ Figure 6(c) ] , thermal stresses due to welding are 
almost zero. The stress distribution along Section B-B is shown in Figure 
6(c) -(2) , Stresses in areas underneath the welding arc are close to zero, 
because molten metal does not support loads. Stresses in areas somewhat 
away from the arc are compressive, because the expansion of these areas is 
restrained by surrounding areas which are heated to lower temperatures. Since 
the temperatures of these areas are quite high and the yield strength of the 
material is low, stresses in these areas are as high as the yield strength of the 
material at corresponding temperatures. The amount of compressive stress 
increases with increasing distance from the weld or with decreasing tempera- 
ture. However, stresses in areas away from the weld are tensile and balance 
with compressive stresses in areas near the weld. In other words, 

fa • dy = 0 (1) 

X 

across Section BB. ** The stress distribution along Section BB is shown in 
Figure 6(c) - (2) , 

Stresses which are distributed along Section C-G are shown in Figure 
6(c) -(3) . Since the weld metal and base metal regions near the weld have 
cooled, they try to shrink, which causes tensile stresses in areas close to the 
weld. As the distance from the weld increases, the stresses first change to 
compressive and then become tensile. 

Figure 6(c) - (4) shows the stress distribution along Section D-D. High 
tensile stresses are produced in areas near the weld, while compressive 
stresses are produced in areas away from the weld. The distribution of 
residual stresses that remain after welding is completed is shown in the figure. 
The stress distribution is quite similar to that shown in Figure 1(a) . 

The cross-hatched area MM' in Figure 6(a) shows the region where 
plastic deformation occurs during the welding thermal cycle. The cross-hatched 


*In a general three-dimensional stress field, six stress components, 

a , cr , cr , r , r , t , exist, 
x y z xy zy zx 

** Equation ( 1) neglects the effect of cr^ and r ^ on the equilibrium 


condition. 
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area near the origin 0 indicates the region where the metal is melted. The 
region outside the cross-hatched area remains elastic during the entire welding 
thermal cycle. 

The above discussion shows that to analyze residual stress and distor- 
tion completely by mathematics, the following steps must be taken: 

1) Analyze heat flow. 

2) Analyze thermal stresses during welding and determine the distri- 
bution of incompatible strains that are produced during the welding 
thermal cycle. 

3) Determine residual stresses and distortion produced by the incom- 
patible strains. 

A number of articles have been published on the subject of heat flow in 
weldments (Section IV) . Although not an easy problem, it can be handled 
analytically. 

The problem of determining the distribution of incompatible strains is 
extremely difficult, however. When the material undergoes plastic deforma- 
tion, the stress-strain relationship is not linear. Furthermore, plastic 
properties of the material change with temperature. 

When the distribution of incompatible strains is known, be it analytically 
or experimentally, the third step, determination of residual stress and distor- 
tion, can be handled analytically. Moriguchi developed a fundamental theory of 
stresses caused by incompatible strains. 24 Masubuchi has applied this theory 
to the study of residual stresses and distortion in weldments. 25 

Because of the difficulty in determining the distribution of incompatible 
strains, no analysis has yet been developed to trace the change of two-dimen- 
sional thermal stresses during welding and to determine distributions of three 
residual-stress components, o^, cr^, and r^. In other words, no analysis has 

been made in which both heat flow and stress fields are treated as two-dimen- 
sional problems. In all studies conducted so far, the problem has been simpli- 
fied in some way. 


b. Analysis of Stresses Caused by Longitudinal Welding in Strips 
and Rectangular Plates 


In 1936, Boulton and Lance-Martin calculated the residual 
stresses produced by depositing weld metals on one or both edges of long, flat 
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plates, 9 feet by 6 inches by % inch. 26 They assumed that the only stress Com- 
ponent that existed was the longitudinal stress (a in Figure 6) and that it was 

X 

constant in the longitudinal direction. Rodgers and Fetcher in 1938 and Rosenthal 
and Zabrs'' in 1940 conducted similar studies. 27 > 28 Since then many investigators 
have studied stresses due to longitudinal welding in strips and rectangular 
plates (including bead welding on plate edges, bead welding along the plate 
center, and butt welding) by making various assumptions concerning stress 
and temperature distribution. 13 ’ 19)29,30 

1) For stress distribution, consider the longitudinal stress only and 

neglect the transverse and shear stresses (cr^ = = 0 in Figure 6) . 

In some cases, it is even assumed that the longitudinal stress is 
constant in the longitudinal direction; or is a function of y only 
in Figure 6. 

2) For temperature distribution, assume that temperature does not 
vary in the longitudinal direction. In some cases, it is even assumed 
that a narrow band in the weld area is heated to a constant tempera- 
ture, while other areas are kept at the initial room temperature. 

Both Tall 31 ’ 32 and Tsuji 33-36 conducted rather comprehensive analyses. 

Tall used the two-dimensional temperature distribution, but neglected trans- 
verse and shear stress (cr = t = 0) . Schematic curves shown in Figure 6 

y xy 

are drawn based on calculations by Tall. 

Tsuji used a temperature-distribution constant in the longitudinal direc- 
tion, as shown in Figure 7. He 
assumed that the central region 
(width 2a) of a rectangular plate 
is suddenly heated to temperature 
T 0 and that heat flows into 
other areas by conduction. He 
handled the plate as a two-dimen- 
sional, elastic-plastic field, 
i. e. , he considered the three 

stress components, cr , a , t . 

x y xy 

Tsuji studied the effects on 
residual stresses of various 
parameters including heat input 
and material properties. 

A recent trend is an in- 
creased use of computers for 
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FIGURE 7. INITIAL TEMPERATURE CON- 
DITION USED BY TSUJI FOR THE 
ANALYSIS OF RESIDUAL STRESSES IN A 
BUTT WELD 
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analysis, especially digital computers. Tall and Tsuji both used computers for 
their calculations. N. N. Prokhorov and N. Nikol. Prokohorov used a digital 
computer for calculating distortion which develops when a bead is deposited on 
the edge of a strip. 37 Makhnenko and Velikowanenko used a computer for calcu- 
lating stresses and distortion in structures made of welded narrow plates. 38 

During the present study, Battelle has developed a FORTRAN program for 
calculating both thermal stresses during welding and residual stresses in a butt 
weld. The analysis is essentially the same as that of Tall. Details of the 
Battelle analysis are described in Section V. 


c. Analysis of Stresses Which Remain After Welding is Completed 


A relatively simple analysis can be made if we are interested in 
stresses which remain after welding is completed and not necessarily in the 
stress change during welding. Based on Moriguchi's fundamental study of 
elasticity and plasticity, Kihara and Masubuchi proposed an approach that can be 
applied to a study of residual stresses in practical weldments. 39 Here, the 
problem is divided into the following two steps: 

1) Determination of distribution of incompatible strains produced by 
welding 

2) Determination of residual stresses due to the incompatible strains. 

The first step can be achieved analytically or experimentally. The 
analytical achievement of the first step is extremely difficult, as described 
earlier. However, it is relatively easy to calculate residual stresses once the 
distribution of incompatible strain is determined. Masubuchi and associates 
have studied various problems related to residual stresses and distortion in 
weldments, which include the following: 

1) Effect of welding procedures on residual stresses and transverse 
shrinkage of simple butt welds 

2) Effect of welding sequence on residual stresses and transverse 
shrinkage of butt welds and circular patch welds. 

Some of the results are described in a paper written by Masubuchi. 25 

The approach that Masubuchi and associates have taken is very similar 
to Prandtl's wing theory. 24 > 1 2 M°> 41 The following is an excerpt that compares 
equations describing distributions of residual stress and shrinkage of a welded 
joint and the phenomena that occur in the velocity distribution around the wing of 
an airplane. 25 
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Theory of Residual Stress 



Stress: a , <j , t 
x y a 

Slit length (weld length) : L 

Transverse Stress on Slit: a . 

yo 


= JE 1/2 1 

^yO 47 r / X - X* 
-L/2 


d[v]\ 
dx yx 


dx' 


where: 

E = Young's modulus, 
[v] = dislocation. 


( 1 ) 


00 

[v] = /, A sin n0 
LJ „ n 
n=l 


(2) 



_E_ 

2L 


Z 

n=l 


nA 


n 


sin nfl 
sin 9 


( 3 ) 


x = — cos Q 

Lt 
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Mean Dislocation: [v] 


1/2 

m = f / [y]dx 
L - 1/2 



sin 0 • d0 



(4) 


The magnitude of mean dislocation is determined only by 
the value A j. 

Strain Energy: U 


where: 


U = 



2 °y0 


[v]dx 




= ^E(l+ 6) A/ 



(5) 


(5’) 


The magnitude of strain energy becomes minimum when 
<5 =0 provided that the length of L and mean dislocation [v] 
remain constant. When <5=0, the form of dislocation [v] 
is elliptical and the distribution of transverse stress a 
is uniform along the slit. ^ 
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Wing Theory 



Velocity: u, v 
Wing breadth: b 

Down Wash: w 


w = 


_1 

47T 


-m x ~ x ' 



x’ 


dx 


i 


where: 

r =» circulation. 


T = 2bV J] A sin n 9 
n 


n=l 


w 


= V Yj nA 


n=l 


sin n8 
n sin 9 


V = velocity at infinity 


x = — cos 9 
z 
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Lift: L 


b/2 

L = pV f (x) dx 
-b/2 

= |pb 2 V 2 Ai 


where: 
p = density. 

The magnitude of life is determined only by the 
value Aj. 

Induced Resistance: D. 


b/2 

D. = p f T(x)w(x) dx 
1 -b/2 



/ 00 \ 
pb 2 V 2 1 Yj nA 2 n 

\ n=1 y 


= — pb 2 V 2 (l + 6)Ai 2 


where: 


<5 = 


1 

A.2 

l 



The magnitude of induced resistance becomes 
minimum when 8 = 0, provided that the wing breadth 
and lift L remain constant. When 6=0, the distribu- 
tion of circulation T is elliptical and the distribution 
of down wash is uniform along the wing. 
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Equation ( 1) in the above excerpt shows the relationship between the 
dislocation* which occurs along the part of the x-axis between (-1/2, 0) and 
(1/2, 0) and the transverse stress produced at that part of the x-axis designated 
as cr . The physical meaning of this stress problem can be interpreted as 


yo 
follows 


When transverse stress cr^, given by Eq. (3) , is applied along both 


sides of the slit of length L located on the x-axis, the relative displacement of 
the sides of the slit is given by Eq. (2) . The value of U given by Eq. (5) can 
indicate the elastic strain energy stored in the plate (per unit plate thickness) 
when transverse stress cr Q is applied along the slit. 


To determine the distribution of residual stress in a weld joint and to 
investigate the influence of welding procedures on stress distribution, it is not 
always necessary to trace the history of nonelastic deformations. The important 
thing in these problems is to determine the distribution of incompatibility. If 
the distribution of incompatibility is determined, residual stresses and distor- 
tions can be calculated theoretically. This approach makes it possible to study 
residual stresses and distortion analytically without becoming involved in complex 
calculations related to plastic behavior of the material. 


3. Buckling Distortion in Weldments in Thin Plates 

An excellent review of buckling distortion in thin welded plates is 
presented in Chapter 11 of a book by Kihara, Watanabe, Masubuchi, and Satoh, 
published by the Society of Naval Architects of Japan. 13 The following discussion 
is taken primarily from that chapter. Although the discussion is concerned 
with carbon steel welded with covered electrodes and with the submerged arc 
process, it is included in this report because no such comprehensive work on 
aluminum is available. Because this is taken from Japanese work, all 
measurements in this section of the report are in the metric system and 
temperatures are in centigrade. 


*The term dislocation used here is not the same as the dislocation used 
in physical metallurgy, although their fundamental natures are closely related 
mathematically. As used here, dislocation means discontinuity in displacement 
such as opening of a crack in a residual-stress field. In this case, displace- 
ments are expressed in multivalued functions. 
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FIGURE 8. APPARATUS FOR 
DEFLECTION MEASUREMENT 



WELDING CURRENT 100 A 
ARC VOLTAGE 25 V 
TRAVEL SPEED 1.0 cm/ sec 
SIZE OF ELECTRODE 3.2 mm 
IN DIAMETER 


FIGURE 9. DEFLECTION DUE TO 
WELDING VERSUS TIME CURVES 


a. Distributing Features in 
the Distortion of Thin 
Plate Due to Welding 


In welding of a thin 

plate , the plate may buckle due to com- 
pressive thermal stress during cooling* 
as shown in Figure 5(c). Fundamental 
investigations of buckling-type distortion 
were undertaken by Watanabe, Satoh, 42,43 
Masubuchi, 44 and Mura. 45 

To obtain characteristics of distor- 
tion in thin plates. Watanabe and 
Satoh 42 ) 43 observed the distortion 
produced when a single-pass, straight- 
bead weld is laid on the center line of a 
steel sheet 700 millimeters long and 1. 6 
millimeters thick. The apparatus used 
to measure deflection of the plate during 
welding and cooling is shown in 
Figure 8. 43 A cantilever of spring-steel 
sheet was set on a rigid station near a 
specimen and the edge was bent to be in 
contact with the specimen. Thus, deflec- 
tion of the specimen was transmitted to 
the cantilever. The spring-steel sheet 
cantilever was made as thin as possible 
to minimize the effects of its rigidity on 
deflection of the specimen. As shown 
in the figure, an electric wire strain 
gauge was attached to the cantilever 
surface. Changes in reading of the 
gauge were measured by a dynamic wire- 
strain indicator. Then curves showing 
the relation between deflection of the 
specimen and time were obtained. The 
results are shown in Figure 9. 43 

As seen in the figure, an abrupt 
change in deflection in the 700-millimeter 
wide specimens occurred about 4 minutes 
after the start of welding. The change 
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WELDING CURRENT 90 A 
TRAVEL SPEED 0.60 cm/sec 
SIZE OF ELECTRODE 2.6 mm 


FIGURE 10. STATES OF DISTORTION 
PRODUCED AFTER BEAD WELD- 
ING IN A FEW SPECIMENS OF 
MILD-STEEL SHEET 



0 100 200 300 400 500 600 700 

DISTANCE FROM THE START 
OF WELD ALONG WELDLINE (mm) 


BEAD WELDING WAS MADE ON THE 
CENTER LINE OF THE PLATE BY 
SUBMERGED ARC WELDING AT 
410 A, 30 V,AND 2.0 cm/sec TRAVEL 
SPEED . 


was so abrupt that it was accompanied 
by a loud noise. Measurement of 
temperature rise at the middle point of 
the welding bead indicated that that 
point had already cooled to 70° to 80° C 
at that instant. In the narrower speci- 
mens, the abrupt changes of deflection 
were not observed but the deflection 
was much increased after welding 
(Figure 9) , 

Figure 10 shows states of distor- 
tion produced after cooling of specimens 
350 millimeters wide, 700 millimeters 
long, and 1. 6 millimeters thick, 43 
Three different states of distortion were 
observed despite the fact that almost 
the same welding conditions were applied. 
In Specimen No. 1, transverse distor- 
tion was concave and longitudinal dis- 
tortion was convex to the side of the 
deposited bead. In Specimen No. 3, this 
relationship was reversed. A combined 
type of distortion was observed in 
Specimen No. 4, i. e. , the states of 
distortion in transverse and longitudinal 
directions reverse themselves in the 
same plate. 

Masubuchi 44 observed longitudinal 
distortions after cooling of submerged 
arc welded long rectangular plates. A 
weld bead was laid on the center line of 
the 720-centimeter long, 2.3-millimeter 
thick plates. Curves in Figure 11 show 
longitudinal deflections produced in a 
plate 40 centimeters wide. 44 Several 
stable states of distortion are observed 
in the same plate. 


FIGURE 11. SEVERAL STABLE 
TYPES OF LONGITUDINAL 
DISTORTION DUE TO WELD- 
ING IN A THIN PLATE 
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FIGURE 12. BRIDGE-TYPE STRAIN 
GAUGE USED FOR MEASURE- 
MENT OF LONGITUDINAL 
STRAIN DUE TO WELDING 


b. Critical Buckling Stress 
When Welding Thin 
Plate 


Fundamental problems 
relating to buckling distortion in 
thin plates due to welding are 

( 1) critical buckling stress of a 
plate which has both a given 
aspect ratio and thickness and 

(2) the amount of deflection 
after buckling. 
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WELDING CURRENT 100 A 

ARC VOLTAGE 25 V 

TRAVEL SPEED 1.0 cm/sec 

SIZE OF ELECTRODE 3.2 mm IN DIAMETER 


FIGURE 13. LONGITUDINAL STRAIN 
DUE TO WELDING VERSUS TIME 
CURVES 


Watanabe and Satoh deter- 
mined experimentally thermal 
stresses produced during weld- 
ing on the specimens shown in 
Figures 7 and 8. In addition to 
the cantilever for measuring 
distortion, they attached a 
bridge-type strain-measuring 
unit on the specimen shown in 
Figure 12. 43 Active and dummy- 
strain gauges were mounted on a 
steel strip 0. 3 millimeter thick 
and 7 millimeters wide, which 
was bent to the shape shown in 
Figure 12. Then the strip was 
attached to the specimen using 
Screws A and B which were 75 
millimeters apart. The bridge- 
type unit was placed in the longitu- 
dinal direction of the specimen near 
the distortion measuring point shown 
in Figure 8 . A thermocouple also 
was attached to the specimen to 
determine temperature change. 


The strain determined by the bridge-type unit included thermal expan- 
sion in addition to strains that were related to thermal stresses; therefore, 
Watanabe and Satoh called the strain determined by the bridge-type unit 
apparent strain. Figure 13 shows changes of the apparent strain in the 
longitudinal direction e during and after welding. 43 Figure 14 shows tempera- 

X 

ture change 0 and resulting free thermal expansion d0 , where 0 1 is the 
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FREE THERMAL 
EXPANSION (ad) 



0 123456789 

TIME (min) 


WELDING CURRENT 100 A 
ARC VOLTAGE 25 V 
TRAVEL SPEED 1.0 cm/sec 
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SIZE OF ELECTRODE 3.2 mm IN DIAMETER 


FIGURE 14. TEMPERATURE RISE AND 
FREE THERMAL EXPANSION VERSUS 
TIME CURVES 


x IQ 4 



AB s a e 

AC = € (LONGITUDINAL STRAIN) 


& -E x BC, 

Cr 3 -3 

= 21 x 10 d (0.160-0.459) x 10 6 

= -6.3 kg/ mm % 

FIGURE 15. AN EXAMPLE OF THERMAL 
STRESS CALCULATION 


coefficient of linear thermal 

expansion. 43 The following 

equation was used to calculate 

thermal stress in the longitudinal 

direction cr : 
x 

a = E U -ad), (2) 

X X 

where E equals Young's modulus. 
In the above calculation, trans- 
verse stress 0 ^ was neglected. 

Figure 15 shows an example 
of the thermal stress calculation 
for a specimen 700 millimeters 
long and 700 millimeters wide. 43 
The calculation was made at the 
point x= 0, y = 1/4B = 175 mm. 
The amount of compressive 
thermal stress at the instant of 
buckling was determined as 

E x BC in Figure 15. In the 
particular case shown in Figure 
15, cr = -6. 3 kg/mm= (9000 

cr » 

psi) . 


In the narrower specimens , 

no abrupt change of deflection, 

as was shown in Figure 9, was 

observed and the magnitudes of 

a were obtained on the 
cr 

assumption that buckling might 
have taken place at the instant 
shown in Figure 16. Values are 
shown in Table I. 

Compressive thermal 
stress is considered to be 
caused by inherent shrinkage due 
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BUCKLE AT THIS INSTANT 


FIGURE 16. CHANGE OF DEFLEC- 
TION IN A NARROW PLATE 
INCREASING CONTINUOUSLY 
DURING COOLING 


to welding. Thus, critical inherent 
shrinkage of buckling was obtained 
upon the assumption that compressive 
thermal stress is equivalent to stress in 
the plate produced by concentrated com- 
pressive forces P = Eh (| j) acting at 
the both ends of the weld line (Figure 17) . 
When the force approaches the critical 

value P the plate buckles. The stress 
cr 

at the measuring point (x - 0, y = B/4) 
due to the concentrated compressive 
forces can be obtained according to 
the theory of elasticity. When l/B = 1,. 
it is given by 


TABLE I. CRITICAL VALUE OF INHERENT SHRINKAGE 

(£ j) FOR BEAD-ON-PLATE AND NUMERICAL 
s 1 cr 

FACTOR k 


L 

(mm) 

B 

(mm) 

cr 

cr 

(kg/mm 2 ) 

P 

cr 

(T) 

(mm) 

k 

700 

700 

6 . 3 

7.45 

0.222 

67 



5. 9 

6.98 

0.208 

63 

700 

350 

9. 9 

5. 55 

0. B5 

25 



8. 8 

4. 94 

0. 147 

22 

700 

230 

14.7 

5.41 

0. 161 

16 



16.2 

5.96 

0. 177 

17.5 


cr 

X 


at x = 0 , 


y 



m 

1+4 (-1) 2 (m7r M + l)e _m7rM , 

m=2, 4, 6 


( 3 ) 


where ju equals L/B. Critical values of P and (|{) (= P/Eh) can be obtained 
(Table I) by substituting a in Table I in the left-hand side of Eq. (3) . 
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P- Eh (£ ) 

WHERE 1 

(£ )■ INHERENT SHRINKAGE 
h « THICKNESS OF PLATE 

FIGURE 17. LONGITUDINAL DIS- 
TORTION AND STRESS IN THE 
PLATE 

the following equation (Figure 18) 43 : 


Watanabe and Satoh have expressed 
the critical inherent shrinkage as 

( .) = IzL krrV 

■ i; cr Eh 12 (1 - v 2 ) B ’ 

where 

h = plate thickness 
v = Poisson’s ratio. 

The term k is a numerical factor 
depending on the aspect ratio l/B of the 
plate. As seen in Table I, the magni- 
tude of k decreases as L/B increases, 
and it is approximately represented by 


k = 


11 + 50 (l) 


(5) 


and is independent of material. 



0 0.2 0.4 0.6 0.8 1.0 


(B/L) 2 . 


Watanabe and Satoh further 
analyzed the relationship between weld- 
ing conditions and the inherent 'shrinkage 
(| . Based on their previous work, 

they proposed the following formula 
for welding carbon steel using ilmenite- 
type covered electrodes 4 millimeters 
(0. 16 in. ) in diameter: 


Ui) 


0. 136 x 10 -6 



> ( 6 ) 


FIGURE 18. NUMERICAL 
FACTOR k 


where 

(| !) = inherent shrinkage (mm) 
I = welding current ( A) 
v - arc travelling speed 
(cm/sec) 

h = plate thickness (cm) . 
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BREADTH OF PLATE (m) 

FIGURE 19. CRITICAL PLATE 
THICKNESS FOR LOW- 
CARBON STEEL FOR VARIOUS 
WIDTHS OF PLATE AND 
ASPECT RATIOS 


It is considered that a plate will 
buckle in the longitudinal direction when 
the magnitude of .(|j) exceeds its criti- 
cal value (I j) obtained by Eq. (4) 
cx* 

and (5) . Figure 19 shows the relation- 
ship among the critical thickness of the 
steel plate, the aspect ratio, and the 
thickness of the plate under the condi- 
tion of 260, where I is welding 

current in amperes and v is welding 
speed in centimeters per second. 43 
Figure 20 shows the relationship among 
minimum half wavelength welding condi- 
tions and the width of the plate on plates 
1. 6 millimeters thick. 

Mura 45 conducted a mathematical 
analysis of buckling- type distortion of a 
long strip due to welding (Figure 21) . 
The width is 2a, the length is L, and 
the weld is made along the center of the 
strip. The following assumption is 
made regarding residual stress 
components: 



~ i i i i i .i _—i 

I 0 0.5 1.0 1.5 2.0 2.5 3.0 

BREADTH OF PLATE (m) 


Longitudinal stress 
Transverse, cr = 0 

y 

Shear, r = 0. 
xy 


-a' /y/> C 

cr" 0</y/< C 


(7) 


Since residual stress must be balanced, 


a 

/ cr dy - 0. (8) 


FIGURE 20. CRITICAL HALF WAVE- 
LENGTH FOR LOW-CARBON 
STEEL FOR VARIOUS WIDTHS OF 
PLATE AND WELDING 
CONDITIONS 


Therefore, 

a' = — - — or” . (9) 

a - c 
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FIGURE 21. WELD SPECIMEN 
USED IN MURA’S ANALYSIS 


If we denote the deformation in 
the z-direction as W, the equation of 
equilibrium is 


D 


9 4 W 

9x 4 


+ 2 


9 4 W 

9x 2 9y 2 


9 4 W 1 

9y 4 y 


= T 


9 2 W 
x 8x 2 


( 10 ) 


where 


D = rigidity of the plate = 


Eh 3 


12 ( 1 - p 2 ) 


T - a h 
x x 

h = plate thickness. 


For boundary conditions, it was assumed that the plate is simply 
supported at both ends, x = 0 and L, and that the plate is free along y = ± A. 
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Mura determined the 
stable wavelengths by 
solving Eq. (11). Curve 
A of Figure 22 shows the 
relationship between non- 
dimensional parameters, 

A. and p , which express 
wavelength and residual 
stress, 


m 2a 


FIGURE 22. RELATIONSHIP BETWEEN WAVE- 
LENGTH AND RESIDUAL STRESS IN A LONG 
WELDED STRIP 


-If 


12 ( 1-V 2 ) <J ! 


where m is the integer. 

The relationship is given as 


(ID 
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From knowing the residual stress cr j, ju is determined. Then A is 
determined. Wavelength values Lj, L 2 , Lg for M = 1, 2, 3, . . . . The 
longest wave Lj is obtained under ju = 1. For a given residual stress 
CTj and for width a, a plate longer than Li is going to buckle. Therefore, Lj is 
considered as the critical size of the plate to cause buckling due to residual 
stress. 


Curve B in Figure 22 shows the Ji-\ relationship when the compressive 
stress 0 - ! is applied across the entire width of the plate. For the same value of 
A, the value of j u for Curve A is considerably greater than that for Curve B. 
This is because the tensile stress in the region near the weld restricts the 
tendency to buckle caused by compressive stress and the compressive stresses 
must be much higher to cause buckling. 


c. Deflection of Plate After Buckling 

Attempts were made to analyze the amount of deflection of a 
plate after buckling. Masubuchi 44 investigated the relation between the wave- 
length of the longitudinal deflection curve and the amplitude in rectangular plates 
720 centimeters long and 2. 3 millimeters thick where the welding bead was 
deposited on the surface by submerged arc welding. Several stable states of 
distortion are observed in this specimen, and generally the amplitude 
increases as the wavelength increases , or as the number of waves decreases 
(Figure 10). 
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FIGURE 23. EFFECT OF WAVE- 
LENGTH IN LONGITUDINAL 
DEFLECTION CURVES ON THE 
AMPLITUDE 


TABLE II. NUMERICAL 
CONSTANTS n AND C 


l/2b 

n 

C (cm) 

72 

2.08 

0.0041 

36 

2.06 

0. 036 

18 . 

2 . 02 

0. 076 


ture A at x = 1/2 is given by 


Figure 23 shows the relation 
between the amplitude A and ratio of the 
wavelength and the_width of the speci- 
men. 43 The term X is the ratio of wave- 
length versus width of the specimen, L 
is length of the specimen, and 2b is the 
width of the specimen. Bead welds 
were made by submerged arc welding. 
Figure 23 is approximately represented 
by 

A = C \ n , (13) 

where C and n are constants. 

The magnitudes of the constants 
obtained from Figure 23 are shown in 
Table II. Accordingly, it is understood 
that the amplitude is almost proportional 
to the square of the wavelength. 

Watanabe and Satoh 43 have 
expressed the deflection after buckling 
as follows: The curve of deflection 
along the half wavelength after buckling 
is expressed by a sine curve such as 

y = y 0 sin — , (i4) 

where L is the half wavelength after 
buckling and y 0 is the maximum deflec- 
tion in the half wave. Maximum curva- 


A - 

L 2 


or 



( 15 ) 
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FIGURE 24. MAXIMUM RADIUS OF 
CURVATURE OF BUCKLED 
PLATE 


estimated by using Eq. (15) and Figure 24. 43 


Hence, y 0 is proportional to the 
square of the half wavelength. 
Equation ( 15) is of the same 
type as the experimental formula 
(13) presented by Masubuchi. 
Referring to curves of distortion 
due to bead welding on a center 
line of rectangular plates 700 
millimeters long and 1. 6 milli- 
meters thick, maximum curva- 
tures in the longitudinal direc- 
tion or the magnitudes of A in 
Eq. (15) are obtained as in 
Figure 24. As the magnitude of 
A decreases the width of the 
plate increases. Deflection 
after buckling of a rectangular 
plate of a given size can be 
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Section IV. COMPUTER ANALYSIS OF THERMAL STRESSES 
DURING WELDING AND RESULTING RESIDUAL STRESSES 


In an attempt to obtain a better understanding of stress buildup and 
distortion resulting from welding, a computer program (written in FORTRAN 
digital-computer language) has been developed to simulate the welding process. 
This program uses present theoretically derived mathematical relationships 
that were already available. Although, realizing that these existing equations 
relied on certain assumptions, Battelle felt that they were sufficiently reliable 
to develop a useful analysis through correlation with good experimental data. 


The computer program consists of three basic sections and several sub- 
programs. The main sections of the program are 


1) The calculation of temperature distribution created in a plate by a 
moving heat source, which simulates making a butt or edge weld 

2) The calculation of the resulting thermal and residual stress 
distributions 

3) The development of tables and graphs relating time, temperature, 
stress distributions, and position on the plate as direct computer 
printer and plotter output. 


Although the basic relationships and equations used in the first two 
sections come primarily from the doctorial dissertation of Lambert Tall at 
Lehigh University, 31 the Battelle study includes the following modifications 
and advancements: 


1) It is fully computerized and mechanized, including printing of 
numerical results and machine plotting of graphs. 

2) The program is designed for any material. 

3) The effects of temperature on material properties, including 
thermal conductivity, yield strength, modulus of elasticity, 
coefficient of thermal expansion, and specific heat, are considered 
to some extent. 


1. Calculation of Heat Flow During Welding 

Heat flow during arc welding has been studied analytically and 
experimentally by many investigators including Rosenthal, 46 Adams, 47 ’ 48 
Masubuchi, 49 Christensen et al. , 50 and Rykalin. 51 Rykalin prepared an extensive 
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review on this subject. Battelle Defense Metals Information Center Report 172 
presents the state of the art of this subject and discusses some shortcomings 
of the present analysis. 52 

This section presents a brief summary of the state of the art of the 
analysis of heat flow during welding as it is related to the computer program 
developed in this study. 


a. Generation and Dissipation of Welding Heat 

(1) Heat Generated by a Welding Arc . The electric power of 
the welding arc H* contributes by far the largest amount of heat to the arc. The 
term H is defined by 

H = VI W (j/sec) , (16) 


where 


V = arc voltage (V) 

I = arc current (A) . 

The thermal equivalent of the electric power is 0. 24 VI calories per second. 

The energy input (or heat input) of the welding arc, which is given in 
terms of joules per inch of weld, is commonly used to express the intensity of 
the welding heat source 


60 x 2 , 

V 


where v is the travel speed of the welding arc in inches per minute. 


(17) 


(2) Dissipation of Welding Heat. Most of the heat generated by 
the welding arc is supplied to the workpiece under the welding arc and then 
dissipates into further portions of the workpiece. The remainder of the heat is 
dissipated into the electrode by thermal conduction and is lost to the surrounding 
atmosphere by radiation. It is customary to express the amount of heat 


*Other heat sources include the heat caused by chemical reactions in the 
arc atmosphere and the heat caused by transformation of metal. 
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supplied to the workpiece Q, calories per second, as a portion of the thermal 
equivalent of the electric power of the welding arc (or arc energy) 0. 24 VI as 


Q - 77 * 0, 24 VI, (18) 

where 17 is called arc efficiency, and Q may be called the net heat or effective 
thermal power of the welding arc. 

Figure 25, presented by Christensen et al . , 50 shows measured values of 
arc efficiency for various welding processes and materials including GTA and 
GMA welding of aluminum . 52 A logarithmic presentation is used so as to cover 
the range of energies from 300 to 12, 000 calories (1 to 50 kVA) . Accordingly, 
the efficiency is given by the position of the curves or bands representing each 
process, while the slope of these bands by definition will remain constant at 
45 degrees. Values of r, range 66 to 85 percent for the GMA process and 21 to 
48 percent for the GTA process for welding of both steel and aluminum. How- 
ever, according to the information obtained at Marshall Space Flight Center, 
the values of 77 for welding aluminum with GTA and GMA processes are about 80 
to 90 percent. 


b. General Discussion on Mathematical Analysis of Heat Flow in 
Weldments 

Figure 26 shows schematically the temperature distribution in a 
plate on whose surface a weld bead is being laid at a speed v. Curves 1 to 6 
represent isothermal curves on the surface, while the dotted curves represent 
isothermal curves on the transverse section ABCD. The coordinate axis is 
0 -xyz; the origin 0 is on the surface underneath the welding arc, the x-axis lies 
in the direction of welding, and the z-axis is placed in the thickness of the 
plate, downward. 

The fundamental expression for heat conduction is given by the Fourier 
heat-flow equation, 51 ’ 52 


9£ _ /s 2 0 , 8 2 9 , 9 2 e\ 

9t y3x 2 9y 2 9z 2 ) 


(19) 


where 

9 -(h)- 0o = temperature change 
( 5 )= temperature 
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HEAT SUPPLIED TO WORKPIECE Q (cal /sec) 
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FIGURE 25. MEASURED VALUES OF ARC EFFICIENCY FOR VARIOUS 
PROCESSES AND MATERIALS 
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TRAVEL SPEED, 



FIGURE 26. TEMPERATURE DISTRIBUTION IN A PLATE WHEN A WELD BEAD 
IS LAID ON THE SURFACE 




9 0 = initial temperature 
t = time 

K = — thermal diffusivity 
c p 

C = specific heat 
p = density 

X = thermal conductivity. 

The mathematical analysis of heat flow in a weldment is essentially a solution 
of Eq. (19) for a given initial condition (initial temperature distribution) and 
a boundary condition (shape and intensity of the heat source, geometry of the 
weldment, etc.). 


Two features of the heat flow during metal-arc welding are 

1) The heat source moves, usually at a constant speed, on or near 
the surface of the workpiece. 

2) The size of the heat source (welding arc) is small compared to the 
size of the workpiece. 

(1) Quasi- Stationary State and Nonstationary State . The 
process of heat flow in arc welding is divided into three stages. 

1) Heat saturation process in which the temperature in the 
field, moving together with the heat source, continues to 
rise 

2) Quasi-sationary state in which the temperature distribution 
is stationary for a coordinate system that moves with the 
heat source 

3) Temperature leveling-off process after the welding arc is 
extinguished. 

In the quasi-stationary state the mathematical analysis is simple, since 
the problem can be treated as a steady-heat- flow problem for a moving coordi- 
nate. The quasi-stationary state takes place at a point close to the weld that is 
continued for a long period. Most of the mathematical analyses conducted so far 
on heat flow in weldments have been on the temperature distribution in the quasi- 
stationary state. 
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Heat flow is in the nonstationary state in areas near the start (Stage 1) 
and the end (Stage 3) of a weld, even when the weld is made over a considerable 
length. When welding is performed over a short length, the quasi- stationary 
state is never reached. Mathematical analysis of heat flow in the nonstationary 
state is much more complex than that in the quasi- stationary state. 


(2) Effect of the Shape of the Heat Source. The local effect 
principle of the heat- conduction theory shows that the pattern of heat distribution 
of a local source has a substantial effect upon the temperature distribution only 
in the region adjacent to the source. Therefore, the temperature distribution 
of the workpiece in a region reasonably far from the welding arc can be calcu- 
lated with sufficient accuracy by schematizing the pattern of the heat-flux 
distribution. In Figure 26, for example, temperature distribution in areas 
outside, say, Curve 4 can be calculated with sufficient accuracy by assuming 
that the heat source is concentrated in an elementary volume (point source) 
located at Point 0. The difference in temperature at a corresponding point on 
the top and the bottom surface becomes less pronounced as the point is located 
farther away from the welding arc (Curves 5 and 6) . In such cases, the 
temperature calculation may be made on the assumption that the heat source is 
concentrated along Line 00' (line source) ; the temperature distribution is 
reduced to a two-dimensional plane problem. 


(3) Physical Constants (Linear and Nonlinear Theories) . 
Physical properties, including thermal diffusivity k , specific heat c, density p, 
and thermal conductivity X, are not really constants but change with temperature, 
and the extent of the change depends upon the metal concerned. 53 Almost all 
of the mathematical analyses conducted so far on heat flow in weldments have 
been made on the linear theory in which thermal properties 54 are assumed to be 
constants. Such solutions give approximate results. When values of thermal 
properties are considered to change with temperature, Eq. (19) becomes non- 
linear and the mathematical analysis becomes extremely complex. Grosh et 
al. 55 analyzed heat flow in weldments assuming that thermal properties change 
linearly with temperature. However, only limited work has been done on the 
nonlinear analysis of heat flow in weldments . 


c. Temperature Distributions in the Quasi- Stationary State 

Mathematical analyses of heat flow in welding in the quasi- 
stationary state have been made by Boulton and Lance-Martin, 56 Rosenthal and 
Scherbes, 57 and many other investigators. 
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Using a coordinate (|, y, z) that is moving at the same speed as the 
welding arc shown in Figure 26, 


| = x - vt . (20) 

The temperature undergoes no change in the moving coordinate system, thus 
Eq. (19) is expressed as 46 > 51 > 52 


9 2 9 _ 9 2 9 9 2 0 _ £00 

9£ 2 9y 2 9z 2 K 

Now 0 is a function of position (£, y, z) only. 


(21) 


Solutions have been obtained for various cases including: 


Case 1, three-dimensional case, semi-infinite plate (for very thick plate) 

Case 2, three-dimensional case, finite thickness 

Case 3, two-dimensional case, infinite plate (for large thin plate) 

Case 4, two-dimensional case, finite size (for thin plate with finite size) . 


In most cases , solutions are obtained for a point or line-heat source with a 
given intensity. In other words, the size of the heat source is considered to 
be infinitely small. Solutions for a finite-sized heat source can be obtained by 
integrating solutions for a point or line-heat source. 


Complexity of a solution varies, depending upon the boundary condition. 
Solutions for Case 1 due to a point source and Case 3 due to a line source are 
rather simple. Solutions for a finite-sized plate due to a finite- sized heat source 
are very complex. 

After the complexity and usefulness of these solutions were evaluated, 
it was decided that the solution for Case 3 (two-dimensional) due to a line 
source be used for the Battelle computer program. The specimens analyzed 
in this study (aluminum plates 0. 25 by 30 by 144 in. ) are considered to be large 
enough that ( 1) temperature is quite uniform in the thickness direction in most 
areas where analyses are made and (2) the size effect of the plate can be 
neglected. More complex solutions may be used in future programs to improve 
the accuracy of calculations . The solution for Case 3, which is used in the 
Battelle program, is described in the following pages. 
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The temperature change caused by a line source (intensity: q = Q/h, 
where h is plate thickness) which moves along the x-axis of an infinite plate 
is given by Eq. (22) , 


0 = JL. e zX ^ 

2nX 2k 5 



( 22 ) 


where 

r = si | 2 + y 2 

K 0 (z) = modified Bessel function of the second kind and zero order. * 


2. Calculation of Thermal and Residual Stresses 

The present state of the art of the analysis of stresses caused by 
longitudinal welding in strips and rectangular plates is described in Section III. 
After usefulness of available solutions for analyzing the problems being 
investigated in this study was evaluated, it was decided to use the analysis by 
Tall. 31 ’ 32 This analysis (Section III) uses the two-dimensional temperature 
distribution but it only considers the longitudinal stress a and neglects trans- 
verse stress a and shear stress t . Compressive stresses in the longitudinal 
x ny 

direction that exist in areas near the moving arc are believed to be the primary 
factor causing lateral metal movement. Consequently, it was felt that Tail's 
analysis would be good enough at least for an initial analysis. It is strongly 
recommended, however, that the analysis be further developed in a future pro- 
gram to include ail three stress components. 


^Detailed descriptions on the modified Bessel functions are given by 
Watson and McLochlan 58 ’ 59 

<X> -zt 

K 0 (z) = f - - • dt . 

_1 n/ t 2 - 1 

When z is large enough, K 0 (z) can be expressed 
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Tail's analysis employs a step-by-step method for calculating thermal 
and residual stresses that result when a plate is arc welded at either the 
center line or edge by a moving electrode, assuming that the bead is deposited 
at a uniform rate. At each step, the plastic deformation is considered and 
equilibrium accounted for. 

First the temperature distribution around the moving arc is calculated 
by using Eq. (22) . Then the stress field is divided into a set of transverse 
strips (Figure 27) . The time interval represented by the strip width must be 
short enough that the temperature and thermal stress for each increment may be 
regarded as being constant. Since the greatest changes in temperature occur 
immediately after the onset of welding, the time increments should be short at 
first and comparatively long for the time approaching infinity. 

The calculation starts on a strip some distance ahead of the welding arc 
where the temperature change is negligible and stresses are purely elastic. 

Time zero is fixed on the strip. For example, in the calculation shown in the 
fourth section of the appendix, the heat source is located at T = 9 seconds. In 
this particular case, the welding speed is 0. 233 inch per second, as shown in 
the third section of the appendix. In other words, the calculation started at a 
strip 2. 1 inches (0.233 x 9) ahead of the arc. 


I 

( 




K 

h 
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FIGURE 27. DIVIDING THE STRESS FIELD INTO TRANSVERSE 
STRIPS FOR CALCULATING THERMAL AND RESIDUAL STRESSES 
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Thermal stresses along the strip are calculated by the following 
equation: 

a - ot E 0 , (23) 

where 

ct = stress 

o' = coefficient of linear thermal expansion 
E = Young's modulus 
0 = change in temperature. 

The relation a = E’0 is derived by equating the equations for expansion; 

thermal expansion of the strip AH ~ h 0 oi 0 , where h 0 is the width of the strip 

is equal to the expansion due to stressing AH = h 0 . The temperature increment 

for the next time interval, or the next strip, is calculated and the thermal stress 
increment corresponding to it determined. This thermal stress increment is 
added to that existing at zero time, and this addition is limited by the yield 
stress of the material at the actual temperature under consideration. At each 
step, the stresses should be in equilibrium; i, e. , the summation of compressive 
and tensile stresses must be zero and the moment of these forces must total 
zero across the total plate, 

f <7 • dy = 0 (24) 

f a * y • dy = 0 . (25) 

■ Since temperature distribution is symmetric with respect to y, Eq. (25) 
is automatically satisfied. In order to satisfy Eq. (24) , a constant stress 
(compressive or tensile) across the plate is added (Figure 28) , 

It is also assumed that the amount of stress at a given point does not 
exceed the yield stress of the material at the temperature of that point. In 
areas near the welding arc, temperatures are high but stresses are low because 
values of the yield stress are low at high temperatures. In fact, liquidus 
areas very close to the arc hold no stress. 

Each increment has its thermal stress added to that already accumulated 
from the preceding steps with the limitations that the addition must satisfy both 
equilibrium conditions and that the maximum yield stress existing at that 
temperature may not be exceeded. The stress distribution at time infinity gives 
residual stress after welding is completed. 
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a>E>9 + CONSTANT = STRESS 


FIGURE 28. CALCULATION OF THERMAL STRESSES IN A 
TRANSVERSE STRIP 


In conducting the analysis, the time increments, or the widths of the 
strip, need to be small enough to give a reasonably continuous thermal-stress 
pattern. For example, the time intervals used in the calculation given in the 
fourth section of the appendix follow: 

a) 1- second interval for time 0 to 50 seconds 

b) 50-second interval for time 50 to 1000 seconds 

c) 5000-second interval for time 1000 to 3600 seconds and infinity. 

3. The Battel le Computer Program and Results of Computation 

The program developed in this study is to calculate temperature dis- 
tribution and resulting thermal and residual stresses due to a moving heat 
source. The program is written in FORTRAN IV computer language for use 
on a Control Data Corporation (CDC) 6400 computer system including a CalComp 
plotter. The FORTRAN language is common in computer programming and can 
be adapted readily to any other computer that has a FORTRAN compiler pro- 
gram. The appendix presents the program and examples of calculation in the 
following order: 
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a) Main program 

b) Subprogram 

c) Listing of input 

d) Sample listing of output data 

e) Examples of graphical computer-output plots. 

a. The Battelle Computer Program 

( 1) Materials and Welding Parameters . To facilitate a broad 
coverage of possible parameter values , the program is written to be quite 
general. Computations can be made readily for any material (if the material 
properties are known) and for any welding parameters (arc voltage, welding 
current, arc efficiency, and arc travel speed) . These values are fed to the 
program as input data, as described in the third section of the appendix. 

A single value for any of the required material properties (density, 
specific heat, thermal conductivity, coefficient of thermal expansion, Young's 
modulus, and maximum yield stress) is all that is required. However, to 
improve the accuracy of the analysis , a table of up to ten values of each property 
as a function of temperature may be entered as data. The temperature range 
used is in equal increments from 0°F to some value near the melting point. A 
given material property value can thus be determined for any temperature by 
means of a special interpolating subroutine (appendix, second section) . 

The program also allows the width of the plate to be subdivided into two 
areas of different increment spacing to allow smaller numerical integration 
areas close to the heat source. In addition, the length of the plate can be 
divided into three separate time spacings of different increments to allow more 
accurate calculations in the area of rapidly changing temperatures. 

The entire set of input values, requiring nine IBM cards, is listed and 
each separate input value is explained at the beginning of the computer program 
description in the appendix. A typical listing of these input values for a run is 
also listed in the appendix. 

After the program has read in and processed all the input values, the 
time and position spacings are calculated along with limiting values and spacings 
to be used on graphical outputs. 
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(2) Computation of Temperature Distribution . The next 
program step is to calculate a temperature value at a particular point by the 
method previously described using parameter values for either a temperature of 
a previous calculation or for ambient temperature ( 0°F temperature difference) . 
The property values at the newly calculated temperature are determined and 
another new temperature is calculated. This process repeats in an interpre- 
tive fashion until the latest temperature calculated is within 0. 1° of the previous 
temperature value. The temperature is calculated in this fashion at each point 
in the plate as specified by the input spacings. These temperature values are 
printed out in table form in the appendix. They are stored for use in other 
calculations and for graphical outputs. 

The actual equation used to calculate each temperature value is of the 

form 


T = K' e a K 0 (b) . (26) 

Since K' is essentially a constant for any single calculation and a and b 
are functions of distance away from the heat source, the time since the heat 
source passed over any given point on a plate will eventually approach a large 
value, cc. As the distance away from the moving heat source increases, the 
parameters a and b approach an infinitely large value. This results in an 

increasingly small value of e multiplied by an increasingly large value K 0 (b) 
to yield a finite value. Since the computer only handles finite numbers (x) in 
the range 10" 247 < x < 10 +347 , it was necessary to incorporate a routine to 

_ a 

approximate, by series summation, the product e K fi (b) for large values of 
a and b. 


( 3) Computation of Stresses . After all the temperatures at a 
cross section have been calculated, the stress distribution is calculated as 
described previously. This involves computing both the change in temperature 
from the previous cross section and an average temperature over the increment 
to determine the values of limiting material properties, namely yield stress 
and modulus of elasticity. This stress distribution, if totaled across the plate, 
must equal zero to fulfill the required equilibrium conditions, Eq. (24) and (25) , 
especially Eq . (24) . 

First, the values of aE0 along the transverse line are calculated. Since 
the temperature change is greater near the center line of the plate and decreases 
away from the center line, values of at Ed will tend to follow the same pattern 
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and hence not total zero across the plate. Thus, a constant stress is added to 
satisfy the equilibrium condition (Figure 28) . If the stress value in any region 
is higher than the yield stress at the temperature of that region, the stress 
value is further adjusted so that the stress does not exceed the yield stress. 

This analysis, however, has an obvious shortcoming regarding the 
effect of plate width on the stress distribution. * Since the equilibrium amount 
is spread out equally over the width, the stress distribution is directly 
dependent upon the plate width. When the width is infinitely large, as an 
extreme case, the stress becomes zero. There is a definite need for improving 
this analysis; the problem is discussed further in Section VII. 


b. Trial Calculations to Improve the Program and to Establish 
Computation Procedures 

In the initial stages of developing the program, different com- 
puter runs were made incorporating a range of values for heat- input energy, 
welding speed, material thickness, type of material, plate widths, and incre- 
ment spacing between points across the plate (y-axis in this analysis) and time 
intervals along the length of the plate (x-axis in this analysis) . In addition, 
modifications in the logic of the program, as well as some analysis changes, were 
made to improve the program's output results. It was determined that the cost 
of drawing a graph on the computer was as low as $0. 75 and seldom over $1. 50 
per complete labelled plot. Thus, a major emphasis was placed on developing 
graphical computer output in a highly descriptive visual form. As can be seen 
from the graphical computer output shown in the fifth section of the appendix, 
a clear understanding of the complexities of the analyses is attained, especially 
in stress-distribution change during a welding cycle. 

To check the validity of the computer program, calculations were made 
on mild steel welds and compared with results obtained by Tall. 31 > 32 

To check the accuracy of the numerical integration part of the analysis, 
tests were run to determine minimum spacing intervals along both the x- and 
y-axes. It was found that spacings in the y-direction as wide as 0. 5 inch made 


*In the present Battelle analysis, a plate with a finite width is used to 
calculate the stress distribution, while the solution for an infinite plate is used 
to calculate the temperature distribution. It was felt that the solution for the 
infinite plate is good enough for the present study. In future studies, a solution 
for a finite-sized plate may be recommended. 
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little difference in the calculations. Spacings in the x-direction of up to 10 
seconds for the first 50 seconds, 100 seconds for 1000 seconds, and extremely 
large thereafter, only slightly affect the final residual-stress calculation. It 
was subsequently discovered that in spite of the fact that the final calculations 
change insignificantly during these intervals, the local thermal- stress distri- 
bution changes significantly and can thus hide the instantaneous effects . Thus a 
lead time ahead of the heat source and a 1-second interval to 50 seconds, a 50- 
second interval on to 1000 seconds, and a 5000-second interval thereafter to 
any desired stopping point were used. 


c. Examples of Calculations 


To study effects of welding parameters on the temperature and 
stress distributions, this report presents results of calculations for bead-on- 
plate welds in 2014-T6 aluminum made under six different conditions (Table III) . 


TABLE III. WELDING PARAMETERS USED FOR SAMPLE 
CALCULATIONS 


Condition 

Number 

Voltage 

(V) 

Current 

(A) 

Efficiency 

Factor 

Speed 
(in. /sec) 

Thickness 
(in. ) 

1* 

10 

254 

0. 8 

0.233 

0. 25 

2 

25 

250 

0.9 

0. 3 

0.25 

3 

25 

250 

0.5 

0. 3 

0.25 

4 

25 

250 

0 . 1 

0. 3 

0.25 

5 

25 

250 

0.5 

0. 6 

0.25 

6 

25 

250 

0.5 

0. 3 

0.5 


^Simulates the conditions used for experiments performed at NASA. 


Condition 1 simulates the conditions used for experiments performed at 
the Marshall Space Flight Center, NASA. For Conditions 2 through 6, changes 
were made in thermal power of the welding arc, arc travel speed, and plate 
thickness. The effective thermal power of the welding arc is changed by varying 
any or all of the values of arc voltage, welding current, and arc efficiency 
[Eq. (18) ] . In this study, a moderately high voltage and amperage for actual 
welding conditions were selected and held constant. Only the efficiency value 


47 




was changed to vary the effective thermal power of the welding arc. It was 
assumed that the weld was made along the longitudinal center line of a plate 8 
inches wide and infinitely long. 


( 1) Temperature and Stress Distributions Across the Plate. 
Figures A-l through A- 10 illustrate the temperature and stress distributions 
across the plate at the instantaneous time given on each pair of figures for 
Condition 1. For this time sequence of figures, the heat source is located at 
the 9-second time position. Numerical values of temperature and stress at 
various times are listed in the appendix. Similar results for Conditions 2 
through 6 were obtained but are not presented in this report. 

Figures A- 1 and A-2 show the temperature and stress distributions, 
respectively, at 7 seconds, or along a cross section slightly ahead of the 
moving heat source. The temperature increase above ambient is 277°F 
at the center line, y = 0. Compressive stresses are produced in regions near 
the center line because the expansions of these areas are restrained by the 
surrounding areas. Tensile stresses exist in areas away from the weld. The 
maximum compressive stress is about 27,000 pounds per square inch, which is 
smaller than the yield stress.* Consequently, the entire cross section is still 
in the elastic condition. 

Figures A-3 and A-4 show the temperature and stress distributions, 
respectively, at 9 seconds, or along the cross section passing the heat source. 
The amount of stress at the center line is almost zero, because the temperature 
is very high. As the distance from the center line increases, the absolute value 
of the compressive stress increases since the temperature of the point decreases. 

At the time increment 15 seconds, the temperature change spreads out- 
ward (Figure A-5) . Stresses in areas near the center line change to tensile 
as these areas start to cool (Figure A-6) . Compressive stresses exist in 
areas between y = 0. 6 and 2 inches. 

At the time increment 1000 seconds, the temperature change is quite 
uniform across the entire section (Figure A-9) . Figure A-10 shows that high 
tensile stresses are produced in areas near the weld while compressive stresses 
exist in areas away from the weld. 


^According to the input data on Card 5, the yield stress of 2014-T6 
alluminum alloy is 62, 000 pounds per square inch at room temperature, and 
it is as high as 50, 000 pounds per square inch at a temperature increase of 
300° F (appendix). 
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The fourth section of the appendix lists temperature and stress values 
for times up to infinity. The stress distribution at infinity represents the 
residual stress after welding is completed. The computer readouts in this 
section show that the stress values do not change significantly after some time 
increment, say 500 seconds. If we take the stress value at the weld center, for 
example, the stress value at infinity is 7, 29 kips per square inch. The stress 
at the center reaches 99 percent of the residual stress, or 7. 22 kips per square 
inch, at 500 seconds, and 95 percent (6.93 psi) at 200 seconds. 

Figures A- 11 and A- 12 show changes of temperature and stress, 
respectively, of the center line (y = 0) , y = 1 inch, and the plate edge (y = 4 in.). 
The temperature at a point on the center line starts to decrease rapidly 
as soon as the arc passes the point, while the temperature at a point on the 
plate edge continues to rise for a long time. Figure A- 12 illustrates that a 
point on the center line is subjected to a complex stress change. First com- 
pressive stresses are produced as the heat source approaches the point. The 
absolute value of the compressive stress first increases as the heat source 
approaches and then decreases as the temperature rises and the yield stress 
decreases. As the heat source passes the point and the temperature begins to 
decrease, the stress changes to tensile and the value increases again. The 
stress at y = 1 inch remains compressive during the entire period, and the 
stress at the plate edge remains tensile. 


(2) Isotherm and Isostress Curves Around the Heat Source . 
Figures A- 13 and A- 14 show curves for equal temperature and equal stress 
around the moving heat source under Condition 1. Figures A- 15 through A-24 
show isotherm and isostress curves around the heat source under Conditions 2 
through 6. Figures A- 15, A- 17, A- 19, A-21, and A- 2 3 show isotherm curves, 
while Figures A- 16, A- 18, A-20, A-22, and A-24 show isostress curves. 

The general pattern of the isostress curves is schematically shown in 
Figure 29. Portions of a plate start to experience significant stress changes 
when the heat source approaches them. There is a high compressive stress 
area shortly ahead of the heat source; however , values of compressive stresses 
are very low in high-temperature areas very close to the heat source. There is 
a narrow tensile-stress zone along the center line behind the heat source, and 
the compressive-stress zone expands outwards from this (Figure 29). 

From the standpoint of metal movement in the thickness direction during 
welding, with which the present study is primarily concerned, the compressive- 
stress zone appears to be of particular importance. If the plate is thin, or the 
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FIGURE 29. TENSILE- AND COMPRESSIVE-STRESS ZONES AROUND 
THE MOVING HEAT SOURCE 


compressive-stress zone is large, the plate may buckle during welding. Dis- 
tortion in the thickness direction also may be caused if the temperature and 
stress distributions are not uniform in the thickness direction. This problem 
will be further discussed in Section VI. 


d. Effects of Welding Parameters 

on Temperature and Stress Distributions 

As described previously, the current analysis has several 
shortcomings, one of which is the way in which the effect of plate width is 
handled. Another shortcoming is that no experimental verification of stress 
values has been obtained. Consequently, results of calculations presented in 
Figures A-l through A-24 should be regarded as showing trends. 
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Figures A- 16, A- 18, and A-20 show isostress stress curves under 
Conditions 2 (high heat energy) , 3 (medium heat energy) , and 4 (low heat 
energy) . By comparing these figures, one could deduce that the following 
phenomena would take place when the heat energy is reduced: 

1) Maximum tensile residual stress in the weld center may not be 
reduced significantly, but the width of the tensile-stress zone 
would be reduced. 

2) Consequently, the amount of compressive stresses and the size 
of the compressive-stress zone could be reduced. 


e. Assumptions and Limitations Inherent in the Present 

Battelle Program 

The presently available information regarding thermal stresses 
during arc welding is very scarce. Consequently, the Battelle program provides 
valuable information on this subject. However, the present program involves 
several assumptions, thus the usefulness of the program is still limited. It 
is hoped that the program can be further developed in the future and its useful- 
ness expanded. The following are assumptions inherent in the present 
program: 

1) The thermal constants and the material properties at 
elevated temperatures are known. (Actually, the data 
available are very limited. ) 

2) The temperature and stress changes in the thickness 
direction can be neglected. (Actually, bending distortion 
due to the stress changes in the thickness direction could 
be very important for causing lateral metal movement. ) 

3) The effects of transverse stress can be neglected. 

(Actually, the weld while melting and cooling causes 
changes not only in the longitudinal direction, but also 
transversely, i. e. , a moment exists due to the transverse 
forces which is neglected in the computation. This moment 
can be increased or decreased if adjoining parts cause a 
nonuniform restraint, for example, at the junction of a 
flange and web. The effect, however, is of secondary 
importance. ) 
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4) The fiber has no further resistance to load when its stress 
reaches the yield-point value. This assumes an idealized 
elastic-plastic stress-strain relationship at all tempera- 
tures. The stress- strain curve at elevated temperatures 

is not of the idealized elastic-plastic type, but rather a non- 
linear curve. The yield point is therefore defined by an off- 
set so that considerable resistance to load can exist past 
the specified yield point. 

5) The heat input is known. (Actually, heat input for this 
investigation was not determined experimentally. The heat 
input has an effect on the thermal and residual stress 
distribution. ) 

In the present analysis, it was also assumed that any transverse 
sections remain plane during the entire process. It was necessary to simplify 
the analysis so that a two-dimensional stress field is replaced by a one- 
dimensional field system. The equilibrium condition also is simplified. 

Instead of using an ordinary equilibrium condition in a two-dimensional stress 
field, simple conditions [Eq, (24) and (25) ] are used. The equilibrium force 
is distributed evenly over the cross section. Consequently, the residual 
stresses could be zero when the plate width becomes infinite. 
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Section V. EXPERIMENTAL INVESTIGATIONS CONDUCTED 
AT THE MARSHALL SPACE FLIGHT CENTER 


Experimental investigations were conducted at the Welding Branch, 
Manufacturing Engineering Laboratory (MEL) of George C, Marshall Space 
Flight Center, NASA, on metal movement during welding of aluminum plates. 
The experiments are described in two reports prepared by R. W. Jackson, 
Welding Development Group of Hayes International Corporation MEL Opera- 
tions. 60,61 The MEL Technical Report WD-70- 6 7 60 describes the results on the 
first test (Panel No. 1) , while the MEL Technical Report WD-79-68 61 
describes the results on the second test (Panel No. 2) . 

This section presents the results from these reports that are important 
to the Battelle study. 


1, Experimental Methods 

Metal movements during welding were measured on two plates of 
2014-T6 aluminum alloy in the horizontal position. The movements were deter- 
mined by 32 dial indicators and recorded on 17 time-sequence photographs. The 
test conditions including welding conditions of the two panels were essentially 
the same. 


a. Test Panels 


Figure 30 shows the test panels and locations of dial indicators. 
The test plates were 144 inches long, 30 inches wide, and 0. 25 inch thick. The 
panels were suspended lengthwise in an upright position between two stanchions. 
The end nearest the weld starting point was securely attached to the stanchion. 
Two cables, two ratchet action puller hoists, "come-alongs, " and two 
dynanometers were used in securing the end near the weld stopping point to the 
other stanchion. This permitted free movement of the panel during welding. 

A tension preload of 3500 inch-pounds was applied to the panel with the come- 
alongs as measured in inch-pounds by the dynanometers (capacity of 0 to 
20, 000 in. -lb) . Tension was maintained throughout the welding cycle. 


b. Welding 

The weld extended for 90 inches along the horizontal center of 
the panel, beginning and ending 27 inches from the ends (Figure 30) . The panel 
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was welded with a bead-on-plate technique using the gas tungsten-arc process. 
The following welding conditions were used: 


1) Arc current 254 A 

2) Arc voltage 10 V 

3) Arc travel speed 14 in. /min (0. 233 in. /sec) . 
The weld was made on the side away from the dial indicators. 


c. Measurement of Metal Movement 


Metal movement was measured with 32 dial indicators mounted 
on a stationary frame in 4 vertical rows of 8 each so that the probes were in 
intimate contact with the surface of the panel opposite the torch. These rows, 
designated A, B, C, and D, intersected the horizontal or weld center line at 
points 6, 30, 54, and 78 inches from the weld starting point (Figure 30) . The 
dial indicators were arranged vertically along these rows at points 1, 2, 4, and 
8 inches on each side of the weld center line. These points were designated 
1 through 8 from the top to the bottom. For example, Point 5 is 1 inch below 
the weld center line. 

Data were collected during tests by time-sequence photographs beginning 
at the time the arc was struck and continuing at intervals until torch shut-off. 
The time interval between frames was about 30 seconds. Another photograph 
was taken approximately 10 minutes after torch shut-off. Each photograph 
simultaneously recorded 32 measurements. The visual reading of the dial 
indicators was made approximately 3 hours after the test with the panel at 
room temperature. 


2. Test Results 

Data obtained by the 32 dial indicators were analyzed to determine 
the metal movement during welding. The following pages describe typical 
results pertinent to the Battelle study. 


a. Metal Movement on Panel No, 2 


Figure 31 shows how points at dial indicators No. 5 (1 in. below 
the weld center line) on Rows A through D moved as welding progressed on 
Panel No. 2. Positive values indicate movement toward the welding torch; 
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FIGURE 31. TYPICAL PROFILES OF METAL MOVEMENT IN THE HORIZONTAL PLATE 




negative values indicate movement away from the welding torch. Four marks, 

A, B, C, and D, in the figure indicate the time when the welding torch passed 
Rows A, B, C, and D, respectively. 

Values of metal movement first increased and then decreased gradually 
as welding progressed. This is shown by curves for Rows A through D. There 
was a definite time lag between a curve for a certain row and that of the next 
row. This indicates that there was a local metal movement near the welding 
torch that traveled with the welding torch. 

The figure also shows that there was a gradual increase in the amount of 
metal movement as the welding progressed from Row A toward Row D. 

Figure 32 shows profiles of vertical axes along Rows A through D at time 
frame 9. At this time the welding torch was located iust beyond Row C. The 
metal movement along Row B was quite significant, while little metal movement 
occurred along Row D. 

These figures indicate two types of metal movement; (1) local metal 
movement around the welding arc and (2) general metal movement of the entire 
panel. 


The local metal movement is insignificant in areas ahead of the welding 
arc. It reaches the maximum some distance behind the arc and then decreases 
gradually. In the particular case shown in Figure 31, the local metal movement 
reduced to an insignificant value about 300 seconds after the arc had passed. 

Or, stated differently, the local movement occurs in areas around 7 inches long. 
It is believed that the local metal movement is associated primarily with the 
thermal stresses around the moving arc. The local movement occurs toward 
the welding torch. This might mean that the plate surface near the torch is 
heated to a higher temperature compared to the other surface; thus the 
expansion of the torch-side surface is greater than the other surface. 

The general movement of the panel increases as welding progresses, or 
as the total heat input to the panel increases. 


b. Metal Movement on Panel No, 1 


Figures 33 and 34 show typical profiles of movement in the hori- 
zontal line and vertical lines, respectively, of Panel No. 1. Figures 33 and 34 
are presented in the same way as Figures 31 and 32. 
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DIAL INDICATOR SITES 



PROFILE OF MOVEMENT ON THE VERTICAL AXIS AT TIME FRAME 9, 
COINCIDENTAL WITH INDICATOR POSITIONS A, B, C, AND D. 

0 ZERO REFERENCE VERTICAL PROFILE BEFORE WELDING. 
A, B, C, AND D VERTICAL PROFILE CHANGE AT TIME FRAME 9. 

+ MOVEMENT AWAY FROM TORCH. 

- MOVEMENT TOWARD TORCH. 


FIGURE 32. TYPICAL PROFILES OF MOVEMENT IN THE VERTICAL AXIS 
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HORIZONTAL PLANE 







Results shown in Figures 33 and 34 are similar to those shown in Figures 
31 and 32, except that results shown in Figure 33 fluctuate considerably com- 
pared with those shown in Figure 31 and the amount of metal movement on 
Panel No. 1 was somewhat less than that on Panel No. 2. 

During the experiment on Panel No. 1, experimental procedures, 
including techniques for measuring metal movement, were not completely 
established. Consequently, results on Panel No. 1 may contain more error 
than those on Panel No. 2, or the amount of metal movement may actually vary 
from one panel to another. Further study is needed to clarify these differences. 

By using the data obtained on Panel No. 1, a series of figures was 
drawn at the Marshall Space Flight Center showing the shape of the panel at 
various instances during the welding cycle. 


Figure 35 is an example of these figures. The x mark shows the loca- 
tion of the welding torch at this instance. The torch is moving from left to 
right. The broken lines show the original shape of the panel and reference 
lines. This figure shows that two types of metal movement exist, (1) the 
local metal movement toward the torch in areas near the torch and (2) the wavy- 
shaped general movement of the entire panel. 
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Section VI. EXPLANATION OF EXPERIMENTAL RESULTS 
ON THE BASIS OF ANALYTICAL STUDIES 


Metal movement, as described in Section V, is a complicated dynamic 
phenomenon, and no analytical information has been developed which is capable 
of quantitatively analyzing it. As a part of the present Battel le study, an attempt 
was made to study the mechanisms of metal movement in a semiquantitative 
manner, based on the present analysis. 

The analysis was made on the following two subjects: 

1) Local metal movement in areas near the welding arc, to determine 
if the local metal movement can be explained on the basis of thermal 
stresses around the moving heat source 

2) General metal movement of the entire panel, to determine if the 
panel may buckle due to residual stresses. 

Areas near the welding arc experience complex stress changes as the arc 
travels. A computer program has been developed in this study to calculate 
temperature and stress distributions around a moving heat source. Figures 
A-l through A- 14, for example, show results for welding parameters of Condi- 
tion 1 in Table III, which simulates the experimental conditions of Panels No. 1 
and 2. * These results, although they neglect stress changes in the thickness 
direction, are useful for studying the general trend. 


The following are two logical possible causes of the local metal 
movement: 

1) The local metal movement by bending. The temperature and stress 
distributions are not uniform in the thickness direction, thus 
bending movement is created. 

2) The local metal movement by buckling. The compressive-stress zone 
near the heat source is large, or the compressive stresses are 
high, thus the panel buckles locally. 


*In the calculation, the arc efficiency was estimated to be 0. 8. The 
plate width for the calculation is 8 inches . 
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1. Possibility of Local Metal Movement Due to Bending 

The possibility of local metal movement due to bending appears to 

be high. 

According to Figure 31, metal movement increased significantly when 
the welding arc approached the point of interest and it took about 300 seconds 
for this movement to decrease. Computer- calculated stress changes around 
the heat source show a similar trend. The calculated stresses start to change 
significantly when the heat source approaches the point of interest, and it takes 
about 300 seconds for the stress distribution to become almost equal to the 
residual stress at time infinity (Figure A- 12) . 

Figure 26 shows that temperature distribution around the welding arc 
is not uniform in the thickness direction, especially in areas close to the arc. 
Consequently, if stress calculations were made using temperature distributions 
shown by Curves 1, 2, 3, ... of Figure 26, the bending moment resulting from 
nonuniform thermal stresses may be large enough to cause local bending of the 
plate around the arc. It is recommended that further analysis be made on this 
subject. 


2. Possibility of Local Metal Movement Due to Buckling 

The analysis of local buckling of a plate due to thermal stresses 
caused by a moving heat source is extremely difficult. Consequently, only a 
simple analysis was made in the present study to see whether or not local metal 
movement due to buckling is a possibility. As it has turned out, the possibility 
of local metal movement due to buckling is remote. 


A compressive-stress zone exists on both sides of the weld line 
(Figure 29) . According to Figure A- 14, the compressive-stress zone extends 
from y = 0. 6 to 3.6 inches, and the stress values are around 7000 pounds per 
square inch. Therefore, an analysis was made to determine if a simply 
supported panel, 3 to 7 inches wide, would buckle under a uniform compressive 
load on the order of 7000 pounds per square inch. 

The critical compressive stress c r of a simply supported panel under 

cr 

uniaxial loading (Figure 36) is 62 


a 

cr 


= K • 


Ett 2 

12 (1-v 2 ) 



(27) 
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FIGURE 36. SIMPLY SUPPORTED m 

RECTANGULAR PLATE UNDER E 

UNIAXIAL COMPRESSIVE v 

STRESSES 


= - + f) 

(27-a) 


JlPQ 

ii 

ca 

(27-b) 


plate thickness 
panel length 
panel width 
number of wave 
Young's modulus 
Poisson's ratio. 


For a narrow panel (/? > 4) , K = 4; 
therefore, 


a = Ett 2 / h_\ 2 
°cr 3 ( 1- v 2 ) ’ ( B J 


(28) 


for 2014- T6 aluminum plate 0. 25 inch thick (E ~ 10 7 psi for temperatures up to 
about 400 °F and v = 0. 3) 


B (in.) 


Vr (pSl) 


3 250, 000 (well over the yield stress) 

7 46,000- 


In other words , extremely high compressive stresses must be present to cause 
the local buckling. Consequently, the possibility of lateral movement due to 
buckling is considered to be remote. 


3. Buckling of The Entire Panel 

The buckling of the entire panel after welding can be analyzed by 
Mura's method described in Section III [Eq. (11) and Figure 22] . 

For a panel where L = 144 inches and 2a = 30 inches, the value of X 
for m = 1 is 9. 6. Using Curve A in Figure 22, the critical p" is about 1. 8. Then 
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the value of (Tj is only 830 pounds per square inch. The compressive residual 
stresses in areas from the weld of an 8-inch wide plate is about 7000 pounds 
per square inch (Figure A- 10). The mean value of the compressive stresses 
of a 30-inch wide plate is considered to be on the order of 1800 pounds per 
square inch, which is much higher than the above critical value of 830 pounds 
per square inch. Consequently, it is believed that the general metal movement 
of the entire panel is caused by buckling due to compressive residual stresses 
that exist in the base plate region away from the weld line. 
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Section VII. SUMMARY AND CONCLUSIONS 


After identifying weld distortion problems in Saturn V structures, a 
literature survey was conducted to obtain useful information on thermal stresses 
during welding and buckling after welding. Using this information, a computer 
program was developed to calculate the temperature distribution and resulting 
thermal stress created by a moving heat source such as exists when making a 
butt or edge weld in a plate. The computer program is a mathematical model 
of a simple welding process as it affects material properties, temperature, and 
stress distributions in a plate. 

The program is written in FORTRAN IV computer language for use on a 
CDC 6400 computer system including a CalComp plotter. Computations can 
readily be made for any material (provided that the material properties are 
known) and for any welding parameters (arc voltage, welding current, arc 
efficiency, and arc travel speed) . The program is so designed that output 
results are printed and figures are plotted. 

To study effects of welding parameters on the temperature and stress 
distributions, calculations were made for bead-on-plate welds in 2014-T6 
aluminum alloy under six different conditions. Computer-plotted results are 
shown in Figures A-l through A-34. 


An attempt also was made to study the mechanisms of the metal move- 
ment observed during bead-on-plate welding of a flat panel in 2014-T6 aluminum 
alloy. The experiments were conducted at the Marshall Space Flight Center. 

The following conclusions have been drawn: 

1) The local metal movement appears to be due to bending moment 
induced by nonuniformity in the thickness direction of thermal 
stresses in areas near the weld. The local metal movement does 
not appear to be due to buckling induced by compressive thermal 
stresses. 

2) The general metal movement of the entire panel is believed to be 
caused by buckling due to compressive residual stresses which 
exist in the base plate region away from the weld. 

It is also recognized that the computer program developed during this 
study could be profitably used in other stress-oriented programs. In addition, 
the success of this program indicates that the general idea of a computer aided 
analytical approach to experimental work has considerable merit and should be 
considered in other areas. 
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Section VIII. RECOMMENDED FUTURE PROGRAMS 


Although the present analysis provides only an indicative picture of 
stress distribution, the development of such a program is in itself an advance- 
ment in the welding technology. It is recommended that further studies be 
conducted to experimentally verify the program, obtain additional experimental 
data needed for the analysis, and improve the accuracy of analysis. 

1) Experimental Verification. First of all, it is important to conduct 
a series of experiments to study the validity of the present 
analysis, which has not been compared with experimental results. 
Most urgently needed is information on thermal stresses during 
welding. 

2) Assessment of Data Needed for Analysis. The analysis requires 
accurate data on properties of metals (thermal conductivity, yield 
stress, specific heat, ... , etc. , at various temperatures) and on 
arc efficiency. However, such data are scarce. It is recommended 
that experimental studies be made to obtain this type of information. 

3) Improvement of Analysis, As described in the last part of Section 
IV, the present analysis involves several assumptions that limit the 
program's usefulness. It is recommended that future studies be 
made to improve the analysis. For example, temperature and 
stress changes in the thickness direction need to be included to 
analyze the metal movement during welding. The effects of local 
transverse stress near, the heat source also need to be considered. 

Figure 37 outlines a possible specific study drawn along these lines to 
further develop the present computer program. 
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FIGURE 37. SUGGESTED STUDY TO FURTHER DEVELOP 
COMPUTER PROGRAM 
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Appendix 

FORTRAN PROGRAM TO CALCULATE TEMPERATURE DISTRIBUTION 
AND RESULTING THERMAL AND RESIDUAL STRESSES 
DUE TO A MOVING HEAT SOURCE WHICH SIMULATES WELDING 
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MAIN PROGRAM 


The FORTRAN program as listed directly by the computer is given. 

At the beginning of the program is a write-up of all the input cards required 
and the parameter associated with each value on an input card. In all, nine 
cards are required per run, but for this listing a repeating loop (in computer 
language a "DO" loop) was added to process up to six runs at a time. 

It is important to note that this listing was run on a CDC 6400 computer 
system, and although FORTRAN is a generally used computer language, this 
program may have to be modified slightly to run on a different type computer 
system. 

The program, as run on the CDC 6400 computer system, requires 
60, 000 octal-core storage positions. Depending on the number of spacings used 
and the number of output graphs to be drawn, it can process an average set of 
data in 30 to 45 seconds of central processor (CP) time and 50 to 120 seconds 
of peripheral processing time. The cost of each graph as drawn by the computer 
ranges from 75 cents to a few dollars depending on the complexity and number 
of data points involved. 

Special subroutines or subprograms written especially to be used with 
this program are listed in "Subprograms. " 
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SAMPLE LISTING OF OUTPUT DATA 


This section of the appendix gives a sample output listing from the 
computer program. Temperature and stress distributions across the plate 
at different time intervals in the process are listed. The temperature values 
listed are a change from ambient conditions and are always positive. The 
stress values are positive for tensile stresses and negative for compressive 
stresses. These values form a table of temperature and stress that a particu- 
lar point in a plate undergoes over a complete welding cycle or the values at a 
specific point at any instant during the welding process. 

The values of temperature and stress start at the center line and con- 
tinue out from it at the increments listed at the start of the program. Although 
it takes more than one line to list the values, each succeeding line is just a 
continuation of the preceding line for values extending out across the plate. 

An infinite time (after the heat source has passed a plane) computation 
of stress is listed and is the resulting residual stress created by the welding 
cycle. 
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LISTING OF INPUT DATA CARDS 


This section presents a list of input-data values as read into the 
program for one specific run and an example of typically used values. 

The parameter associated with each value is listed at the beginning of 
the FORTRAN program at the beginning of this appendix. 

Cards 1 and 3 contain limits, spacings , and parameters for graphing. 

Card 2 contains, in order, the following parameters: voltage; 
amperage; efficiency factor; welding speed; plate thickness; temperature- spacing 
increment for tables of material properties; melting temperature; a factor to 
increase the yield-stress limit in the plastic zone; and the maximum time to be 
plotted on the contour graphs. 


Cards 4 through 9 contain material-property variation with temperature 
data from 0°F to some incremental value near the melting temperature. Each 
of these cards contain ten values spaced at equal temperature intervals. For 
example, for aluminum the melting temperature is 950° to 1100°F. Therefore, 
the temperature increment may be 100 a F, ranging from 0 s to 900° F. 

The properties on cards 4 through 9 in order are 

Card 4: thermal conductivity (W/in. -°F) 

Card 5: maximum yield stress (ksi) 

Card 6: modulus of elasticity (ksi x 10 3 ) 

Card 7: coefficient of thermal expansion (l/° F) 

Card 8: specific heat (W-sec/lb-°F) 

Card 9: density (lb/in. 3 ). 

The card number is punched in the card in column 1, the next parameter 
in the next 7 columns , and one parameter in every 8 columns to the end of the 
card (card number plus 10 parameter values) . 

An example of these input data cards is shown on the following page. 
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SUBPROGRAMS 


The three subroutines listed in this section are used frequently by the 
main program and may be called on from any point in the main program. 

The first subprogram (UART) determines the value of a given material 
property (N) at a given temperature value (Temp) . The routine selects two 
values in the input table, between which the given temperature lies, and 
performs a linear interpolation. This method was selected over many others, 
including fitting an equation to the data, because of the discontinuities of actual 
material properties as a function of temperature. These discontinuities 
generally result from phase changes and change in crystal structure during 
temperature change. 

The second subprogram (BESMKO) calculates the value of a Bessel 
function and is used in the temperature calculation. 

The third program simply labels the computer drawn graphs with a 
legend of descriptive values for the specific calculation. 

Other mathematical and plotter routines were used but these were 
called directly from the computer’s library. Since these programs are 
industry standards or company private, they are not listed. 
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STRESS 44,15 45.43 27.10 -.23 -13.77 -13.09 -12.31 -11.41 -10.40 -9.30 -8.12 -6.89 




o 

rH 

4 

n 

rv 

tn 

ro 

P~ 

40 

O 

rM 

IV 

o 

tv 

4 

CO 

tv 

4 

rv 

© 

c 

40 

r-f 

ro 

04 

o 

• 

• 

• 

» 

• 

• 

• 

* 

» 

• 

• 

• 

• 

ro 

O 

H 

c 

€2 

c 

Cr 

o 

fx 

o 

40 

O' 

in 

fx 

rH 

fx 

rH 

f- 


X 

0—* 

40 

rH 

40 

« 

x 

0-4 

• 

•H 

1 

»H 

1 

r-1 

t 

H 

1 

—9 


rH 


ro 

o 

CO tv 

© a* 

40 4 

tv O' 

in 

40 

oo co 

to © 

h- o 

ro 4 

r-t rH 

r> © 

O' 

fx 

O' 

co O' 

o tv 

M fx 

ro in 

rH 

m 

>o CO 

o* tv 

o o 

O' co 

in rH 

n i w 

O' 

• 

• 

• • 

• • 

• • 

• » 

• 

• 

• • 

• p 

• * 

• • 

• P 

• * 

• 

O 

tv 

«-< ro 

O' 4 

-* ro 

f- m 

rH 

m 

in x 

o ro 

4 co 

© tv 

(V O' 

© (V 

© 

cb 

40 

w-4 

c- x 

rH 

|4- 40 

•HI 


r- 40 


fx X 

r~i 

|X 40 

0-4 

rx 

rH 


« 

r- 4 

1 

0-4 

f 


>-« 

1 

0-4 

• 

•9 

1 

•9 


CO 0—4 

in n 

O © 

O' X 

p- rv 

in © 

0-4 sC- 

rn 

© ro 

ao © 

co to 

in ►- 

tv 

rn r> 

tv co 

4 40 

O' 40 

<3- CO 

O' ro 

40 O' 

fx rH 

CO o 

X co 

o o 

4 x 

co 

• • 

• • 

:■ * 

• ♦ 

• • 

• • 

• • 

* * • 

• • 

• • 

• • 

• • 

■ 

fx rH 

tv rv 

in tv 

•h rv 

ro to 

h rv 

*H ^ 

•h rv 

O' x 

rH *h 

co fx 

0-4 rH 

-o 

ao p- 

rH 

CO h- 

0-4 

CO I- 

rH 

cc r- 

0-4 

r- r>< 

•H 

fx h- 

*»9 

p- 

rH 

» 

0-4 

1 

rH 

1 

0-4 

1 

tH 

t 

•9 

t 

rH 


40 O 

© CO 

in p- 

tv to 

tv o 

in co 

in fx 

co in 

in © 

tv (V 

tv fx 

Kh 

4- 

•h rv 

© v0 

O' co 

r-- in 

CD 4 

X rH 

K ^ 

ro © 

Fx nf 

ro fx 

CO (V 

O 40 

O' 

ro © 

CO 0—4 

O rH 

tv r^ 

cr tv 

0J rH 

x to 

tv 0—4 

•t ^ 

rv 

tv in 

tv 

© 

O' co 

0-4 

O' CO 

rH 

CC' CO 

w—4 

a- co 

H 

CO 00 

rH 

00 CO 

rH 

00 

rH 

« 

0—4 

1 

w-4 

1 

•H 

l 

0-4 

f 

rH 

r 

*H 


© CO 

IT CC 

O' or* 

4 in 

CO © 

in tv 

9* 4 

in o' 

O' rH 

co © 

© rH 

0-4 O' 

O' 

O 4* 

X 4# 

in 4 

to to 

CV 4 

tv © 

o tv 

O' © 

CC © 

CC 4 

C0 FX 

4© ^ 

Fx 

03 O' 

n 

in © 

ro 

m 0-* 

to 

rH OJ 

tv 1 

00 ro 

tv »’ 

x rn 

tv r 

4 

o* co 

0-4 

O' O' 

•-4 

O' O' 

0-4 

O' O' 

rH 

CC O' 

rH 

CO O' 

rH 

CD 

rH 

1 

rH 

1 

rH 

t 

0-* 

1 

rH 

9 

rH 

• 

rH 


tv to 

© X' 

x. ro 

4 X 

© ro 

in 

X 

4 FX 

x in 

tn 4 

in fO 

in in 

4* O' 

rv 

® © 

rH |X 

tv QC 

(V CO 

a* in 

ro 


9* rH 

x tv 

r-* fx 

fx m 

O' tv 

© in 

CD 

rH ©4 

« 1 

O' O' 

» 

x © 


rH 

4 rH 

rH 

IV rH 

rH 

oTv 

rH rH 

fx 

© O' 


O' o 


a © 


t 

O' © 

1 

O' © 

# 

cc © 

I 

CC 

tv 


rH 


0-4 rH 



*H rH 


H r-4 


rH oH 


•H 


o 

tv 

in © 

tv O' 

© CO 

tn co 

tv in 

O' rH 

in rH 

rH CO 

m x 

tv © 

in © 

tv 

in 

X 

tv O 

O' H 

X o 

m x 

fx m 

00 rH 

© 

in x 

rH X) 

tv CO 

V s fx 

© 

4 

cc 

r- tv- 

*-9 O 

tv- 

O' CP 

tV 

X cr. 

to rv 

© 

co tv- 

tv © 

co tv 

o 

© 

tv 

© 

rH 

tv i 

© o 

tv rH 

tv 1 

O' © 

rH rH 

rv 9 

O' rH 

rH rH 

tv 1 

O' rH 
•H r—4 

tv # 

O* rH 
rH rH 

tv • 

O' 

•H 


fx 

tn 

O' rx 

ro fx 

in rH 

© r~> 

CC 

in 

fx ro 

rH CC 

-9- 00 

tv © 

© CD 

tn tv 

4«f 

tv 

ro 

in rv 

in x 

O' ro 

O' O' 

© 

in 

fO rH 

uo 

O' tv 

in cc 

x tn 

CO CD 

tn 

1 • 

:* 

• • 

• • 

• • 

• • 

* 

• 

• • 

» • 

• « 

• * 

• # 

• • 

• 

X 

® 

in tn 

ro co 

in tn 

© cc 

X 

ro 

ao o' 

x rn 

in o 

x tn 

tn O' 

x tn 

r-t 

© 

rH 

-9- • 

© rH 

9f 1 

© rH 

4* 

l 

O' r9 

49- 1 

O' rH 

<r i 

O' rH 

i 

O' 

tv 

rH 


fV rH 


fV rH 



rH rH 


rH rH 


rH 0-4 


•rH 


to in in in in 

o o o o o 


tv 

ro 

tv 

Z 

X CO 

cc tn 

Z 

H CO 

r-l 

in 

z 

x tn 

in 

in 

Z 

rH 4 

X 4 

Z 

in rH 

»x m 

at 

rH 

ro it 

o 

© » « 

x m 

O 

*4 rH 

cr 

■|X 

o 

CO rH 

r-H 

rx 

o 

4 © 

CV O' 

O 

© O' 

in O' 

* • 

.• 

• * 

u 

• • 

• • 

u 

• * 

• 

• 

o 

• * 

P 

r 

o 

• • 

p p 

o 

p p 

p # 

X 

CO 

49 >9 

u? 

>9 CO 

49 4 

IU 

•H CC 

•4 

4 

UJ 

co a 

Xi 

4 

U-’ 

X CO 

*n 4 

UJ 

4 ix 

in 4 

© 

tv 

4 1 

o. 

© tv 

4 1 

in 

C- (V 

4 

i 

in 

o ty 

4 

8 

in 

O' OJ 

4 I 

in 

O' CV 

4 1 




rv 0—4 



tV 









rH rH 



rn 



4 



in 



X 



|x 



tn 

4 


tr> 

4 


in 

4 


in 

4 


in 

4 


in 

in 


a 

tn 


a 

tn 


n 

in 


C 

in 


a 

in 

UJ 


5 

Ul 


> 

UJ 


> 

UJ 


>_ 

Ul 


X 

UJ 

cr 


IU 

nr 


UJ 

a 


10 

tt 


UJ 

cr 


UJ 

QC 

H 


F* 

F— 



*-■ 


F— 

*~ 


»- 

F— 


1- 

fx 

in 

ii 


m 

1) 


in 

li 


in 

tt 


tn 

II 


tn 


u* 



U- 



U' 



UJ 



UJ 







>. 



"S 



X 



X 










i 



M 



►H 




*- 



F- 



»- 



9— 



F- 




94 


TIME = 48 SECONDS 





o 

o 

O' 

>4 

OJ 

4- 

O' 

to 

%o 

CO 

r-t 

0J 

fVJ 

CD 


O' 

.ft- 

?- 


cc 

OJ 

OJ 

<x> 

ft- 

in 

>0 

o 

in 

in 


• 

• 

• 

• 

• 

• 

« 

• 

• 

V 

• 

• 

• 

. • 


O' 

to 

O' 

OJ 

a 

fH 

CO 

o 

ft- 

CO 

K 

o 


to 


1 

xD 

I 

vO 

rt 

j 

0J 

rt 

1 

o 

Ht 

l 

CO 

C 

CO 

i 

i 

in 

4 X© 

© 4 

-4 co 

in t- 

x© r-t 

»© r- 

>4 OJ 

1- o 

O T- 

in *4 

O' O' 

*H CD 

O' CO 

x© 

Ht 

ft- to 

in nt 

in M 

o o 

in o 

«-» -4 

in in 

CO o 

O' © 

oitn 

to f- 

in « 

in nt 

* 

m 

• m 

• • 

• • 

• 

* • 

• • 

.♦ • 

• * 

• • 

• • 

• • 

:» » 

• • 

• 

© 

© OI 

O' rt 

© oj 

CO OJ 

© nt 

to 4 

co to 

Hi O 

ft- in 

O' 4 

o- in 

o o 

r- xO 

to 

> 

rt 

x© »- 

«-» 

>o 0- 


OJ CO 

t l 

O GO 

i « 

COf- 

i » 

CD vD 

1 ft 

ft- 


ft- 

n in 

>c in 

OJ 01 

4 ft- 

H*t to 

4 O 

nO © 

O' 

4 to 

x© CO 

OJ OJ 

O x© 

IV O' 

CO OJ 

ft- o 

• 

• • 

• • 

» • 

• • 

• • 

• * 

• ■ 


*-» nt 

4 CD 

Ht rt 

to o> 

r— t 

-4 CD 

CO 4 

ft- 

r-t 

K N 

•H 

ft- ft- 

r-t 

0J co 

» t 


1 

Ht 

1 

r-t 

1 

» 




m o 

OJ «-• 

o o 

x© o 

1- O' 

co in 

4 .0* • 

o 

o to 

Hi co 

CO 0J 

to in 

ft- o 

OJ co 

O' 4 

• 

• • 

• • 

• • 

• • 

• » 

• • 

* • 

%o 

OJ 

o> o 

r-c 

ft- r- 

-« » 

JO *Ht 

CO 4 

CO 

rt 

0- CO 

Ht 

ft- cc 

r-t 

OJ O' 

1 * 


9 

*-t 

» 

Ht 

1 

- Ht 



x© 

in o* 

4 4 

© in 

in f- 

in oj 

>D O' 

hi to 

«o 

in n 

CO O' 

tO CD 

O' 4 

0J — 

4 to 

r-t O' 

• 

• • 

• • 

• • 

• ® 

• • 

• • 

* • 

4 

OJ ft 

0J 4 

OJ ft 

o in 

OJ Ht 

r- m 

O' 4 

O' 

r-t 

CO O' 

H 

CD O' 

1 

OJO' 

ft 1 


ft 

r-t 

ft 

r-l 

« 

Ht 



o 

X© 

ft— Ht 

CC Ht 

00 x© 

ro xr 

O' OJ 

0J CO 

ft- 

HI CO 

ft- rt 

to O' 

ft- 4 

4 rt 

ro cc 

cd ro 

• 

• • 

» • 

• * 

• • 

• • 

• • 

• • 

0J 

rt H 

in to 

Ht Ht 

ro ro 

rt OJ 

CC CO 

4 in 

o 

Ht 

1 

co o 

Ht rt 

1 

CD © 

Ht Ht 

ft 

Oj O' 

Ht 

a 


co in 

OJ x© 

ro co 

ft- 

O' 

ro ro 

4 ro 

CO 

ft- OJ 

o ro 

co o 

r- 

O' 

Ox co 

x© ro 

OJ 

• • 

• .« 

• :» 

• 

• 

• » 

• • 

• 

OJ 0.1 

co m 

O' xO 

ft- 

in 

© o 

1- x© 

4 

© CO 

» i 

00 ft- 

1 

> 

CO ft- 

i ft* 

1— 


xO OJ 

OJ 4 

oj r- 

ro co 

O' ro 

O' ® 

t© 

in 4 

r-t ,© 

ro in 

CD r-t 

OJ O' 

X© 4 

m 

* * 

• * 

• • 

• • 

• • 

• • 

• 

ro 4 

co in 

©r- 

ft- xO 

Ht Ht 

ft- xD 

4 

© cc 

ft ft 

o ft- 

1 1 

00 ft- 

* * 

ft- 


© Ht 

OJ rt 

ro O' 

o 

ft- 

CO CO 

4 Ht 

O' 

oj in 

oj a 

ft- O' 

O' 

ro 

in o' 

ft- XO 

ft- 

• • 

• a 

.* • 

• 

• 

• • 

• • 

• 

4 xO 

cc u> 

© CD 

f- 

X© 

Ht OJ 

ft- X© 

4 

O 00 

ft ft 

O'?- 

< 

f 

co r- 

ft ft 

ft- 


© OJ 

ft- f- 

4 in 

H 4 

Ht O' 

00 4 

X© 

ft- in 

r- *-t 

© ro 

Ht in 

00 O' 

Oi ft- 


m • 

• • 

• • 

• • 

• • 

• • 

• 

4 00 

in xp 

Ht O 

X© X© 

Ht IO 

X© X© 

4 

© CC 

» 

O' cc 

ft 

a ft- 

» 

ft- 



fO 4 

O' cc 

Ht ?— 

ro x© 

ft- 

© 

X© ft- 

4 rt 

X© ro 

m oj 

ft- 4 

O' Ht 

x© 4 

ft- X© 

CD 

O 

m O' 

CO OJ 

® O' 

00 4 


OJ 

O H-t 

4 ® 

© 4 

4 4 

0J so 

ft- r- 

O' O' 

Ox CO 

© 

• 

• • 

• • 

• • 

• • 

• 

• 

• • 

• • 

* • 

• • 

• • 

• • 

• » 

• • 

• 

Ht 

« OJ 

ft- Ht 

CC OJ 

IP Ht 

ct 

ro 

O' OI 

O IT. 

in © 

ro x© 

#— r-t 

fO xD 

Ht 4 

ro x© 

in 

rt 

OJ « 

ep rt 

OI 1 

«0 Ht 

OJ 

t 

OJ © 

ro i 

© O' 

ro « 

a co 

ro i 

CO ft- 

to ft 

ft- 


4 

Ht fO 

f- x© 

© to 

x© 4 

r- to 

x© © 

in to 

ft- 4 

x© in 

© X© 

to ft- 

r- 4 

Ht co 

X© 

4 

OX O 

Ht 4 

OJ © 

© 4 

CV OJ 

4 4 

O' OJ 

oj o« 

O' x© 

4 00 

to cc 

© 00 

in O' 

H» 

• 

• • 

• • 

• * 

.* » 

• » 

• * 

9 m 

• • 

• • 

• • 

• • 

• • 

• 9 

• 

O' 

X© 4 

O' O' 

?- 4 

1- O' 

ft- 4 

O' in 

© X© 

in oj 

Ht X© 

rt OJ 

OJ X© 

oj in 

OJ xO 

in 

rt 

Ht 

4 1 

a> Ht 

rt nt 

4 | 

CO rt 

Ht rt 

4 t 

OI © 

rt rt 

in i 

© O' 

Ht 

in f 

O' CD 

in • 

C0 ft- 

in f 

r- 




in 



in 



in 



in 



in 



in 





o 



o 



o 



o 



o 



a 



in 

x© rt 

2 

O' x© 

4 CO 

z 

ft- 4 

(V x© 

z 

O' O' 

to to 

z 

in x© 

in ft- 

z 

in nt 

OJ OJ 

z 

Ht® 

Ht O' 

r- 

x© — — t 

O 

in in 

O' © 

o 

4 tO 

© tv 

o 

tn 4 

ft- X© 

o 

to o 

I- 00 

o 

4 © 

Ht o 

o 

Ht X© 

to © 

• 

• • 

u 

• • 

• • 

a 

• • 

• • 

o 

• 9 

» • 

o 

• • 

9 9 

o 

• • 

• • 

o 

• • 

• • 

?- 

in in 

LU 

O' ft- 

in in 

UJ 

ft- ft- 

x© in 

UJ 

O' X 

O' vC 

UJ 

in ro 

© X© 

UJ 

Ht Sj 

rt f»- 

UJ 

OJ X© 

Ht f— 

OJ 

4 ft 

in 

co tv 

4 1 

m 

CC OJ 

4 i 

in 

CV © 

4 9 

Of 

© O' 

in « 

w 

O' CO 

in i 

in 

CO ft- 

in « 




Ht ft 



Ht rt 



rt r*l 



»t 










O' 



© 



o 



© 



© 



© 




V) 

4 


in 

in 


in 

© 


in 

in 


tn 

o 


in 

in 


in 


V) 


Q 

in 


a 

in 

Ht 

a 

in 

rt 

a 

tn 

OJ 

a 

in 

OJ 

a 

in 


bj 


X 

LU 


X 

UJ 


3fc 

UJ 


X 

UJ 


X 

UJ 


X 

UJ 


QC 


UJ 

a 


UJ 

cr 


u. 

a 


UJ 

cr 


UJ 

oc 


UJ 

cr 


ft- 



p- 


ft- 

>* 


»- 

*- 


ft- 

ft— 


ft— 

»- 


ft- 

ft— 


V) 

n 


tn 

a 


in 

II 


in 

II 


in 

u 


in 

ii 


in 



w 



UJ 



UJ 



u- 



UJ 



UJ 





X 



>: 






T 



X 



X 





rt 






rt 



M 



rt 



rt 





ft- 



>- 



ft- 



ft- 



ft- 



ft- 




95 


TIME * 300 SECONOS 



9-1 

to 

ro 


4 

CV 

4 

to 

in 


in 

00 

*0 

00 

CVJ 

o 

CV 

c 

cu 

or 

CVJ 

c 

CVi 

4 

CVJ 

■4 

cv 

31 

• 

• 

• 

♦ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

I*- 

4 

a- 

ro 

t*- 

O' 

t- 

■ h- 

' h- 

4 

hr 

CVJ 

r- 

© 

1 

4 

l 

4 

t 

in 

t 

in 

a 

in 

f 

in 

1 

tn 


O' 

00 

O' 

4 

O' 

4 

O' 

r-* 

CO 

in 

00 

O' 

00 


CVJ 

ro 

cu 

4 

cv 

•-* 

CVJ 

ro 

cv 

00 

cv 

4 

cv 

K 

• 

m 

• 

• 

• 

* 

• 

* 

9 

» 

m 

• 

• 

* 


h- 

f'- 

to 

n- 

o 

r** 

t* 

h- 

4 

h- 

cv 

a- 

o 

» 


• 

4 

i 

4 

i 

in 

1 

in 

« 

in 

1 

in 


in 

o 

4 

O 

CO 

• n 

cv 

4 

(V 

K 

*-» 

CD 

*-# 

♦ 

ro 

CO 

ro 

t- 

to 

4 

to 

in 

ro 

© 

ro 

CO 

m 

»o 

• 

• 

■m 

• 

• 

• 

9 

• 

• 

• 

• 

• 

• 

• 

r- 

r- 

h- 

ro 

a- 

o 

h- 

a- 

a- 

in 


cv 

K- 

© 

a 

4 

1 

4 

t 

4 

I 

in 

1 

in 

rj- 

in 

1 

in 


9* 

r- 

O' 

© 

a- 

© 

4 

O' 

in 

a- 

* 

4 

ro 

o 

4 

H 

ro 

rl 

m 

N 

ro 

a- 

ro 

cv 

fry 

o 

co 


• 

• 

• 

9 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

t- 

CO 

h- 


a- 

© 

a- 


r- 

in 

a- 

ro 

r- 

w* 

♦ 

4 

I 

© 

f 

4 

t 

in 

t 

in 

t 

in 

t 

in 


r- 


to 

00 

H 

to 

O' 

O' 

a- 

in 

4 

cv 

tn 

4 

4 

in 

4 

ro 

4 

o« 

ro 

O' 

ro 

4 

to 

cv 

ro 

CVJ 

* 

• 

9 

• 

• 

• 

• 

• 

.9 

• 

• 

• 

• 

• 

a- 

CO 

h- 

4 

h- 

© 

a- 

a- 

h- 

in 

r- 

ro 

a*- 

r-4 

i 

'O 

1 

© 

1 

4 

9 

in 

1 

in 

• 

in 

i 

in 


tn 

rv ro 

-» rv 

h- O' 

ro O' 

4 r-t 

4 4 

cv 

© 

a- in 

© 

4 

«-t © 

co N 

4 4 

r- 4 

© 

in 

in 4 

CO o 

4 tn 

>0 O' 

4 h- 

«-*to 

4 

CO 

f-* © 

4 

a? 

V) o 

ro o' 

to «-* 

ro O' 

r> 

• 

• • 

• • 

• • 

• • 

9 9 

9 9 

• 

• 

• • 

• 

• 

• * 

9 9 

• • 

• • 

• 

in 

a- 'O 

CO CV 

a* 4 

4 co 

h- sO 

»-* 4 

a- 

4 

00 4 

a- 

4 

in tv 

a- 4 

co o 

a- 4 

r>4 

4 

i a 

4 4 

a a 

4 in 

1 1 

'em 

t 

a 

in in 

« 

i 

in m 

i « 

tn in 

« i 

in 



© in 

co in 

r-» O' 

in © 

f- O' 

n co 

<?“ vO 

ro o> 

cv tv 

^ O' 

© r- 

<0 ^ 

O' © 



in in 

©>- 

in vo 

C0 vO 


ro cc 

^ 00 

ro ^ 

O' 

r- ro 

V O' 

in 

ro © 

♦ 

* 

• • 

• 9 

• • 

• • 

• * 

• • 

• * 

• • 

• • 

• • 

• * 

• ♦ 

• » 

• 

»o 


O' cv 

h~ <x> 

O' 

r- sC 



CO <J- 

r- nT 

in cvj 

n- o 

ro © 

r- tw 

n 


» i 

v> %0 

V 1 

so in 

l I 

>o in 

1 I 

in in 

t 9 

in in 

l t 

in tn 

» > 

tn 



O sD 

o in 

ro co 

ro co 

O' NO 

co r- 

«X> ro 

voty 

co co 

sD sD 


rH 

co * 

©in 

vD 

ro 

V3 vO 

ro ^ 

in r- 

© «-* 

-a- co 

■vf ro 

O' 

O' 

St O' 

co r- 


O 

IT CO 

sf © 

in 

• 

• • 

• • 

:• • 

* * 

• * 

* • 

• • 

• • 

• • 

• • 

• 

* 

• • 

• • 

• 

h- 

K 'O 

0* ro 


in © 

fw sO 

•-« h-‘ 

r- vO 

CD ^ 

r- 'C 

in cv 

n- 

r- 

ro© 



V> 

I 1 

O 'O 

1 t 

nO nO 

1 I 

o in 

t f 

in in 

* a 

min 

i 

l 

tn m 

i • 

in 


ro 


CO 

00 vD 

CC 4 

CO 4 

© 4 

O' 

O' 

I- ro 

r- 

ro 

>c cv 

in no 

*> v£> 

ro O' 

ro 


>0 r- 

r-» 

in cc 

^ r- 

IT O' 

OCT 

4 

O' 

in m 

4 

© 

O' •-* 

4 © 

'C ri 

4 © 

V? 

• 

• • 

• * 

• • 

9 9 

• 9 ■ 

• » 

• 

• 

• * 

9 

• 

• • 

• • 

• • 

• * 

• 

CD 

r- m 

O' 4“ 

r- © 

in © 

h- sD 

*-» N- 

r- 

'O 

co in 

h- 

n- 

in ro 

f- h- 

CO r-4 

r- h- 

H 

<€> 

i « 

OO 

i r 

%C> nD 

1 t 

no tn 

f 

a 

in in 

t 

i 

in in 

1 » 

in m 

i i 

m 


rK 

O 4 

tv in 

4 > in 

© 4> 

•-< © 

«--* O' 

in in 

4 »-« 

co or 

rv in 

© © 

com 

cv cv 

O' 

O' 

ro cc 

O h- 

4 O' 

to cv 

in © 

r- rv 

in © 

4 r- 

in © 

© 4 

4 r-H 

r- 4 

4 

4 

« 

• • 

■ 9 9 

• • 

• 9 

• • 

• • 

• • 

9 9 

• :• 

• 9 

• • 

• • 

• • 

• 

CD 

*> 'O 

O 4 

4> 4> 

in r-< 

4 K 

CO 

4 r- 

co in 

4 h- 

4 ro 

4 

fO rn 

4 K 

•M 


t 

4 4» 

t 

4) 4 

t 

4 if) 

t 

in in 

• 

in tn 

i 

in in 

t 

in 


in 

4 

60 4 
O O' 

CV 4 
r- ro 

4 CVJ 
*-• o 

CO 4 
ro r- 

♦ r— 
rv © 

CD *-t 

in © 
rv #-< 

cvco 

t- © 

co to 
cv 

co r- 
© r- 

© in 
ro *~t 

co ro 
r- n- 

CV 4 
ro -h 

4 

O' 

4 

4 4 
to | 

o in 
4 4 

4 r- 
co t 

IT — 
4 4 

4 

ro a 

r-s CO 
4 IT 

4 r- 

ro a 

cc < 
in in 

4 n- 
co a 

4 T 0 
in in 

4 f- 
ro a 

ro •-< 
in tn 

4 r- 

ro » 

m 


in 

ro in 


CO r-« 

»*• O 


ro tv 

4 ro 


cv © 

© in 


4 CV 

ro 

h- 


h r- 

in co 


© o 

r- © 


4 

ro 

4 O 


h- O' 

h- H 


4 rv 

r- h 


00 r~* 

00 *-< 


r- 4 

cc 

9-1 


9~f © 

IX> *-< 


CD © 

CD CV 


K 

• 

• • 


• • 

• • 


• • 

* 9 


9 # 

• * 


• • 

• 

9 


• • 

9 9 


• o 

» 9 


• 

o 

CV K 


O' in 

rv f- 


in rv 

tv r- 



(V f- 


00 4 

cv 

h- 


4 4 

rv r- 


ro cv 

(V 


r-l 

n- 

in a 


4 4 

in a 


4 4 

in i 


4 in 

in I 


in in 

in 

9 


in in 

in a 


in in 

in a 


in 



in 



CO 



in 



in 




tn 



in 



m 




o 



o 



o 



o 




o 



o 



D 


© 

ro to 

z 

© ro 

© 4 


in 4 

4 CO 

2 

ro n- 

o © 

z 

4 

m 

9-1 

z 

rv in 

in rv 

z 

9-4 m 

h- rn 

Z 


© 

4 — 

o 

CD 4 

in *-< 

o 

4 4 

in ►-« 

o 

CO 4 

4 CV 

o 

r- r- 

4 

tv 

o 

*-> ro 

4 (V 

o 

co rv 

4 CV 

o 

N 

• 

• • 

u 

• o 

• • 

o 

• • 

• « 

o 

.« • 

• • 

o 

» • 

• 


u 

• • 

• » 

o 

• • 

• • 

o 

* 


N 

UJ 

O' vD 

r— * K 

U! 

ifi rv 

^ n- 

UJ 

r-10' 

-< r- 

UJ 

a 4 

*— i 


IU 

4 4 

9-1 S 

UJ 

ro cv 

*-t N 

UJ 

■ 9-1 


in i 

in 

4 4 

If. 1 

in 

4 4 

if. t 

in 

4 in 

in a 

If* 

in ir. 

in 

a 

in 

in m 

in a 

in 

in in 

in a 

in 

in 



© 



© 



© 



© 



© 



o 



o 


in 

in 


in 

o 


in 

in 


in 

© 


in 

in 


tn 

o 


tn 

in 


tn 

ro 

ft 

in 

4 

a 

in 

4 

ti 

in 

m 

a 

in 

in 

0L 

tn 

4 

a 

tn 

4 

a 

UJ 


3 

UJ 


i 

UJ 


> 

UJ 


3 

UJ 


3 

UJ 


3 

UJ 


3 

a 


UJ 

cr 


LJ 

cr 


UJ 

□ 


UJ 

a 


UJ 

DC 


UJ 

ir 


UJ 

f- 


1— 

»— 


K- 

h 





J- 

■a- 



h- 


h- 



►*» 

in 

u 


in 

II 


if. 

11 


in 

1! 


in 

ti 


tn 

n 


in 

it 



UJ 



U.' 



U> 



U* 



U! 



UJ 



UJ 



31 



s: 



V. 



3. 



X 



3 



3 



►-< 









H 



M 



♦— « 



►-t 



t- 






h- 



h 



1- 



»- 






96 



a£> 

p 


4 

\C 

fO 


in 

p 

P 

P 

;■■ o* 

P 

p 

CM 

CO 

ru 

r 

cu 

O' 

cm 

*> 

CM 

4 

CM 

co 

IM 

co 

• 

• 

• 

» 

• 

• 

• 

* 

• 

« 

• 

• . 

• 

• 

P 

cc 

p 

p 

P 

in 

p 

4 

p 

co 

P 

CM 

P 


t 

4 

t 

4 

t 

4 

• 

4 

t 

4 

• 

4 

9 

4 


CO 

4 

CO 

O' 

CO 

CO 

CO 

CO 

CO 

O' 

CD 

© 

CO 

CD 

CM 

o 

CM 

4 

OJ 

o 

CM 

p 

CM 

in 

CM 

ID 

CM 

4 

• 

• 

• 

* 

• 

• 

• 

• 

• 


• 

• 

• 

• 

P 

O' 

P 

p 

p 

sD 


4 

P 

co 

P 

CM 

P- 


• 

* 

« 

4 

t 

4 

1 

4 

1 

■c 

1 

* 

1 

4 


f>4 

o 

o 

ro 

O 

o 

© 

o 

o 

© 

© 


O 

p 

CO' 

CM 

ro 

vO 

CO 

r\j 

m 

O' 

co 

p 

CO 

in 

CO 

in 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

» 

p 

O' 

p 

P 

p 

>0 

p 

4 

p 

ro 

p 

CM 

p* 

H 

*v 

4 

l 

4 

t 

4 

l 

4 

i 

4 

1 

4 


4 

co 

4 

CM 

a£> 

CM 

CM 

rl 

o 

***• 

© 

w~* 

O' 

9-+ 


co 

co 

ro 

P 

ro 

co 

co 

O 

CO 

00 

CO 

*) 

CO 

>o 

• 

* 

• 

♦ 

■* 

• 

• 

• 

• 

• 

• 

* 

• 

• 

p 

O' 

p 

P 

p 

aC 

p 

in 

p 

CO 

p 

CM 

p . 


i 

4 

• 


» 

4 

i 

* 

9 

4 

9 

't 

t 

4 


* 

P 

4 

p 

CO 

CO 

CO 

o 

CM 

O' 

CM 

P 

CM 

CO 

CO 


CO 

CD 

CO 

4 

CO 

*— * 

CO 

00 

CO 

P 

CO 

P 

• 

• 

• 

• 

• 

.» 

• 

m 

♦ 

» 

• 

• 

• 

• 

p 

O' 

p 

P 

p 

'C 

p 

in 

p 

CO 

p 

CM 

P 

H 

1 

4 

• 

4 

1 

4 

1 

4 

I 


9 

4 

» 

* 


in 

NO O' 
CO O' 

00 -* 
in o 

in cm 
n o 

CO 4 
O' NO 

in in 

CO O 

*H O' 

in co 

■4 P~ 

ro o 

O' 

CM 

^ co 
CO o 

p- 

O' *H 

CO o 
CO 

* CO 

00 r-» 

CO rl 
CO IH 

o 

CD 

CO 

4 

P- -o 

ft 9 

O' P 
4 4 

p p 

i t 

p in 

4 -t 

77 

O 4 

4 4 

P' p* 

i t. 

in co 

p- p- 

» K 

CO CM 

p“ p> 
1 9 

CM ^ 
-t + 

9- p- 

t I 

«♦ 

CM 

00 

00 CO 
CO o 

P O' 
>0 OJ 

P 

co o 

sO O 
O O' 

sT CO 
CO o 

O CM 

o -n 

in o 
m — < 

vO in 
cm <*- 

tn,-« 
. co 

* o 

O CO 

>4- CM 
CO •-« 

o <n 
O' CO 

CO •-» 

CD 

CO 

4 

p p 

1 1 

O' P 

•4- 4- 

P P 

1 1 

cr in 

4 4 

p p 

I 1 

'C 

*(9 'S- 

p- p- 

t s 

in ro 

p- p- 

t i 

CV. 1 

p- P- 

9 » 

CV 

+ * 

p- 
9 1 

*— « 

•# 


CM 

rH 

O' P- 
no 

sO 

p- m 

CO o* 
CO c 

f-ft 

P- »*« 
CO 

sD ^ 
NO 00 

vt> CM 
fO 

cm in 
no 

in y* 
CO w 

O' in 
© in 

in in 
co •-* 

'O ro- 
O' in 

xP *> 
CO w< 

•H 

O' 

<P 

p^f- 
< « 

O' P» 
<p 

p- p- 

9 « 

CO M> 

p- r- 

» i 

yO # 
sf sf 

P- P- 
1 1 

in co 
o- ^ 

p- p- 

9 9 

+ CM 
* sp 

p* p- 
* » 

CM r-f 

'P -t 

P* 
ft ft 

•-I 

^p 


o 

*-4 

CM 

r-t 

© CM 

© P- 

O' 

^4- 

cm in 

cc 

in 

p- in 

CO y£> 

y& CO 

P- 

p- 

© © 

P- CO 

in 

'4- 

•P — « 

«C 

CC 

-<P 

CV CO 

r' 

*-* 

p- © 

*o 


CO o 

CO ^-< 

ri P- 

CO 

»— t 

© p- 

CO H 

o* 

O' 

P- P- 

O' 

p- 

P“ K 

CO >o 

p- 

P^ 

yO m 

p- 


in ro 

p^ 

^ CM 

P- 

p- 

CO w 

p- p- 

H 

-4- 

* * 

^4- 

«4- 

1 • 

* * 

l 

t 

-4- 

9 

t 

4 ’P 

« t 

-4- 4 

c 

1 

«♦ *4 

t • 



© 

CO ^ 

cd in 

in aO 

in o 

O P- 

P- IT 

yC CC 

OJ 

CO 

P- O 

P- © 

or 

O' 

co in 

«0 © 

CO 

yC 

y£f 

co o 

vO *-♦ 

cm in 

« -« 

P* CM 

•c. 

*4 

© 

aO h 

r-« O 

aC 

r-4 

o CO 

y£> CM 

O' 

O' 

*> P- 

O' cc 

yO P> 

co 'O 

sD P- 

%o~in 

OP 

in 

<4 

aO P- 

*4 CM 

aD 

p- 

co •-> 

*> P- 

H 

4- 

« 

« ^4 

» 

<4 * 

1 

«4 sP 

t 

<4 

^4 

1 

-4 sp 


i 

-4 <4 

ft 



CM 

CO P- 

CM CC 

in oc 

CO O' 

aO O' 

O' 

<4 

K © 

4- © 

CD w 

o *c 

CC 


vO © 

O' CM 

© 

O' 

CO »H 

O' CM 

(O —• 

CM P- 

CO r-« 


-4 

CO CM 

■4 ru 

CO CM 

CM O 

m 

CM 

© © 

CO CM 

© 

o 

< p- 

O CD 

■4 P- 

a at 

•4 p- 

At 

itj 

4 P 

in * 

<C p- 

W CO 

<4 

P* 

CO CM 

+ P* 

CM 

«4 

CO 9 

4 <4 

CO 1 

* *4 

n i 


*4 

co 1 

•4 4 

co t 

4 -4 

CO 

1 

*4 < 

CO » 

4- 


in 

o © 


<4 

cr 

© *-• 


© o 

r- CM 


r- CM 

fvi M 


o a£> 

CO CO 


cm r ~y 

4 

co 


P 4 

in 4 


»H 


CO OJ 


O' 

<4 

O CM 


CO O' 

O' CM 


CC aO 

O' CM 


4 ro 

O' CM 


CM CM 

O' 

CM 


© #H 

O' CM 


© 

• 

• * 


• 

• 

• ' • 


• • 

• • 


• • 

• • 


• * 

• • 


• • 

• 

• 


• » 

• • 


• 

o 

CM P“ 


o 

CO 

CM P- 


00 aO 

CM P> 


aO in 

CM P~ 


in 4 

CM P- 


4 CO 

CM 

P- 


CO CM 

CM P- 


CM 

in 

in » 



•4 

in t 


4 4 

in t 


4 4 

in i 


4 4 

in » 


4 4 

in 

1 


4 4 

in t . 


4 



in 




in 



in 



in 



in 




in 



in 




o 




o 



o 



o 



o 




o 



a 


p 

O' co 


in 

O' 

o 4 

2r 

HP 

•H -4 

2D 

CM CO 

cm in 

2D 

sD HA 

co in 

z 

cm in 

4 

in 

? 

CD P- 

in in 

2 

CM 

co 

aO cm 

o 

O' 

43 

9- CM 

o 

CO -» 

t- CM 

O 

CC P~ 

P- CM 

O 

4 in 

P- OJ 

o 

CM CO 

p- 

CM 

o 

O CM 

P“ CM 

© 

© 

• 

:• * 

u 

>• 

• 

• • 

o 

• • 

• • 

L> 

• • 

* • 

U 

• • 

• • 

o 

• • 

* 

• 

u 

• • 

• • 

U 

• 

o 

•-< p~ 

u. 

O' 

tc 

*—• P- 

lu* 

cc t- 

r-4 P- 

1C 

*c- in 

r- 

UJ 

in 4 

•H P- 

IM 

4* ro 


p 

UJ 

CO CM 

H P 

UJ 

CM 

in 

in i 

v> 

4 

4 

If) 1 

in 

4 4 

IT: 1 

cn 

4 4 

ir» * 

in 

4 4 

IT. 1 

in 

4 4 

m 

i 

m 

4 4 

in i 

in 

4 



© 



© 



© 



© 



© 



© 


in 

© 

in 

o 


in 

in 


V) 

© 


in 

in 


in 

o 


in 

in 


o 

in 

p 

a 

in 

p 

a 

in 

a 

a 

Wi 

co 

n 

in 

O' 

Ql 

in 

O' 

0L 

in 

© a 

UJ 


X 

UJ 


£ 

u.» 


a 

UJ 


a. 

tu 


X 

u» 


X 

UJ 

H X 

cc 


UJ 

nr 


UJ 

o: 


UJ 

Cl 


UJ 

cr 


UJ 

cr 


UJ 

a 

UJ 

p- 


p 

p 


p 

p • 


p 

p 


p 

p 


p 

p 


p 

p 

p 

<n 

H 


in 

n 


in 

n 


in 

it 


in 

n 


in 

« 


in 

n 


U.' 



U.' 



u 



u; 



U 1 



U.' 



U’ 


x 



>. 



y 



>. 



X 



>; 



?E 


M 



»*4 



*-« 



M 



M 



M 



M 


P 



P 



*- 



P 



P 



P 



p 




















97 




f- 


in 

r- 


fV 


* 



• 



• 


f- 


k 

9- 


1 


0-4 

1 


CO 


ft 

1^ 


CXJ 


* 

rS 


• 


• 

• 


t- 


r- 

f- 


• 


t—* 

t 


O' 


ft 

9- 


rvj 


«t 

rS 


• 


» 

• 


h- 


9- 



» 


r-4 

1 


*-s 


b- 

K 


n 


«t 

r-4 


• 


• 

• 


b- 


9- 

r- 


i 


»H 

t 


<V 


CO 

9- 


ro 


st 

H 


• 


9 

• 




b- 



« 


*->< 

1 

*n 

(V(V 


cc ft 

9- vO 

(V 

m 


si ro 


» 

» 9 


• • 

• - • 

o 

t~b- 


9- 9- 

9- 9- 

s* 

1 1 


0-4 r-4 

1 9 

<v 

mm 


O' ft' 

a? 'ft 

st 

CO r-i 


«t ro 

H »-* 

• 

• • 


• • 

• • 

o 

t- k. 


9- 9- 

9- 9- 

st 

I i 


•— « r*4 

1 9 

00 

* f- 


O' O' 

00 ft 

in 

m *-« 


■t m 

r~i r-l 

• 

• • 


* • 

• • 

o 



9- r- 

9- 9- 

st 

i i 


0-4 H 

i l 

st 

ft O' 


O' o 

O' ft 

»-* 

r>-« 


<r st 

m-4 r-S 

• 

• • 


• • 

• • 

o 

h- 9- 


r- 9- 

t- r- 

«• 

9 t 


*""• *•"« 

i i 

«0 

a ^ 


O' •-* 

'ft 9> 

CO 

ftrv 


*1 -t 

CD **s 

• 

• • 


• • 

• • 

o 

ft N- 


h* 9- 

'ft 9- 

-♦ 

• 


«— » *-s 

« 

CV 

O' tv 


© (V 

f- 9* 

o 

ro tv 


in * 

in ^ 

• 

• • 


• • 

• * 

w" 

<t k 


9- r» 

st 9- 

< 

to e 


r-* r-S 

co t 

m 

in st 


© ro 

r> 9- 

•-* 

O' tv 


in st 

H rS 

• 

• • 


» • 

• • 

H 

MK 


9- 9- 

ro 9- 

♦ 

tn i 


rS 

in i 



(ft 



f** 

m%o 

o 

z 

© sf 

ro 9- 

(V 

I s - (VI 

o 

in -t 

O' »-s 

-0 

• • 

o 

• • 

• • 

iH 

0~ib- 

UJ 

9- r~ 

r*- 


tn • 

(ft 

rH :*-» 

in i 



o 




in 

© 


v? 


tn 

o 

CL 

in 


UJ 

ft 

X 

UJ 


a: 


UJ 

a: 


t- 


f- 

b- 


v> 

II 


(ft 


iii 

X 


98 


sC 

xft 

V 

(V 

• 

• 

w 

fs- 


1 


'ft 


O' 

ru 

• 

• 

(V 

r- 

r— < 

t 


sft 

sft 

O' 

<v 

• 

• 

(V 

n- 

H 

i 


h- 

'ft 

O* 

rv 

• 

• 

(V 

r- 

p— < 

i 


t- 

'ft 

o 

rv 

■■• 

• 

(V 

h- 


l 

r- tv 

'C in 

O' O' 

tv OJ 

• • 

• • 

tv rv 

r- f- 

r « 

S 1 

^ n 

vr m 

O' O' 

cv rv 

* * 

« • 

cv rv 

t- f- 

r-* r-« 

i i 

r- m 

'ft in 

O' o 

rv tv 

• * 

!» • 

rv rv 


rS rS 

i r 

t'- -d* 

cc in 

O' O' 

rv rv 

• • 

» • 

(V tv 

f- f- 

•—* rS 

i i 

1" St 

cc in 

O' O' 

(V 

• • 

• • 

CV (V 

'ft r- 

r-* •"•s 

l 

cc in 

o in 

O' O' 

eg 

• o 

• • 

tv tv 

s* h- 

r*S rH 

m i 



cc in 

St 'ft 



o o 

© (V 



• * 

9 • 



(V (V 

m h- 



r-s H 

in i 


(ft 



(ft 

a 



D 

z 

eo in 

st o 

Z 

o 

o o 

CO tv 

o 

o 

• • 

• .• 

u 

UJ 

cv cu 

*-» K 

Ul 

Oi 

*~t r-S 

in t 

(ft 


© 



© 

o 


(ft 

© 

© 

a. 

(ft 

© 

r-t 

z 

UJ 

ft) 

, »-s 

UJ 

a 

«-s 


9- 

»- 


H 


(ft 

It 

UJ 



UJ 

X 



•X 

►-» 



l-S 

»- 



»- 


r- 

r- 


t' 

(V< 


• 

• 

* 

c 

r> 

c 


l 


h* 

r- 

tv 

r- 

(V 

st 

• 

• 

■■■ » 

o 

n- 

O' 

r~t 

i 



f- 

t~ 


rv 

r- 

(V 


st 

• 

» 


• 

© 

t- 


O' 


1 



r- 



tv 

r- 

(V 


stf 

• 

* 


• 

© 



O' 


1 



cc 

r~ 


(V 

r- 



st 

• 

• 


• 

© 



O' 

*-s 

i 



cc in. 

r- r- 


(V © 

N h- 

(v rv 


St ^ 

* • 

* • 


* • 

© O 

r- f" 


O' O' 

rS r— < 

r i 



CC tft 

r- r- 


rv © 

1^- f- 

rv rv 


st >t 

• • 

• • 


• • 

C 1 * o 

t- V- 


0 O' 


i • 



CO V> 

r- 


rv 

h“ T^* 

(V rv 


st st 

» • 

» • 


• • 

© o 

h- t- 


O' O' 

*-S rS 

• 1 



CC ft) 

O' > 


tv r-f 

r- h- 

CV OJ 


< st 

» m 

m m 


* • 

© © 

f- 0- 


O' O' 

rS iH 

1 1 



CC ft) 

r- 


(V r-* 

r- r- 

t- (V 


st st 

• • 

• • 


• • 

© o 

'ft r- 


o o 

»-S 

• 



CO ft) 



CV «~! 

h- r^ 

St (V 


st st 

• • 

• • 


• • 

© © 

st t- 


o o 

rS H 

m t 



00 tv 

ro t- 


tv 

s r- 

© rv 


•t st 

• • 

• • 


» • 

© © 

m 


O' O' 

r-S •-» 

in i 





m 




o 


cct- 

m rv 

z 

tv *-» 

t- K 

Cft CV 

o 

st st 

• a 

• * 

(J 

• • 

© © 

f-' t- 

UJ 

O' O' 

»— i .••*< 

in » 

(ft 




© 



tn 

© 


a 

m 

© 

Q. 

£ 

UJ 

<rsr 

z 

UJ 

or 

rv 

UJ 

9- 

»~ 


I- 


ift 

II 




u< 




X 




vs 




1- 



r- 

h- 

a 

(V 

■st 

rv 

• 

• 

• 

I- 

CO 

r- 

1 


i 


k 

r*- 

CO 

rv 

st 

cv 

• 

• 

• 

tv 

00 

t- 

* 


1 

K 

r»- 

CO 

rv 

st 

cv 

• 

• 

• 

s- 

a? 

r- 

f 


9 

t- 

fs. 

(0 

cv 

st 

tv 

9 

» 

• 

t - 

CO 

f- 

• 


1 


r- 


r- 

CD 

tv 


st 

(V 

• 


• 

• 

r- 


cc 

K 

i 




r- 


f" ft) 

CC CO 

rv oj 


st st 

rv rv 

• • 


• » 

• • 

|S> f- 


00 CO 

h- r- 

1 9 



1 1 

DO fv 


f- ft 

co so 

cv rv 


st St 

rv rv 

• • 


• * 

* * 

r- r- 


CC a 

S 

i l 



1 9 

co r*- 


t*. ft) 

CO CO 

ru rv 


st -t 

rv cv 

• 9 


• • 

« • 

r- iv 


CO CO 

t- r- 

9 » 



l l 

O' r- 


fw ft 

O' co 

rv rv 


St st 

tv tv 

• • 


» • 

• • 

r- h- 


00 ro 

»w r- 

9 1 



9 9 

ft? 1"- 


cc I- 

ft CO 

r- rv 


•t st 

h- tv 

• • 


■ 9 9 

• • 

ft) fs. 


00 CD 

ft t- 

1 



9 

r- r~ 


CD 1*- 

N CO 

st rv 


* st 

st rv 

• • 


• 9 

• • 

<t 


CC 00 

st 

ro i 



tn t 

ro tv 


cc r- 

m co 

© tv 


•t St 

c rv 

• • 


• » 

« • 

ro t-> 


co a> 

ro t- 

in t 



in t 


(ft 




O 



<o r- 

z 

CO f- 

ro co 

CD (V 

o 

St st 

co rv 

• 9 

o 

9 m 

• • 

«-< JV 

Ul 

co CC 

0-1 f*s- 

in • 

(ft 


m t 


© 



m 

o 


to 

in 

© 

a 

to 

UJ 

ft) 

X 

UJ 

a 

rv 

UJ 

a 

t~ 


t* 

t- 

4ft 

u 


v> 


UJ 




z 








»- 





ft 

O' 


H 


r- 

1 *“l 


rv 


» 

• 


• 


t- 

f- 


•r^ 



» 




ft 

O' 


■0-4 


f^s 

0-4 


rv 


• 

• 


.• 


r» 

t- 





« 




ft> 

O' 


rS 


r- 

•H 


tv 


• 

• 


• 


t- 

b- 


f*s 



• 




ft> 

O' 


0-4 


h- 

0-4 


(V 


• 

• 


• 


t- 

r- 





t 




h> 

O' 


H 


r- 

0-4 


(V 


• 

• 


• 


r- 

b~ 


fs- 



1 




fw ft) 

O' CO 


r-S 


N-f- 

0-4 w—4 


tv 


• • 

9 • 


• 


t- t- 

b~ b- 


t"- 



1 1 




rw ft 

&> CO 


»4 


fs. 

rS r—i 


rv 


• • 

9 » 


• 


f- r- 

t- r~ 


t*. 



» i 




r- ft 

O' ft 


r*S 


t- r- 

H »— < 


(V 


• 9 

• • 


• 



h- 


9- 



9 1 




t“ft 

© ft 


•-4 


t- h- 

CV r-« 


(V 


• • 

• * 


• 



fs- fs- 


f- 



« 1 




ft 

in cc 


#-l 


h- r** 

CO *s 


(V 


9 9 

• • 


• 


t-l^ 

ft h- 





l 




r- ft 

ft O' 


.40* 



in -s 


rv 


• • 

• 9 


• 



st r*. 


t- 



(O « 




i- ft 

tv O' 


0-4 



IPS «-l 


M 


• • 

• • 


• 



ro fs- 


9- 



in t 



<ft 



tn 


D 



Q 


z 

r- ft 

tV O' 

z 

0-4 

o 


O' *-« 

o 

tv 

o 

• • 

• • 

o 

• 

UJ 

t-t- 

r-S h- 

UJ 

b- 

(ft 


in « 

m 


o 



o 


© 


in 

o 


© 

a 

V) 

o 

(L 

rS 

z 

UJ 

ft 

X 

tn 

UJ 

or 

co 

Ul 


I- 

8- 


9- 

ft 


tf) 

H 


U) 



UJ 


X 



X 


»-t 



»-« 


h- 








r* 


CP 


cv? 


rv 


• 


• 


p 


p 


i 


a 


p 


CP 


CM 


CVJ 




* 


P 


p 


a 


a 


p 


CP 


fU 


OJ 


« 


• 


p 


P 


a 


a 


p 


CP 


CVI 


CVI 


♦ 


• 


p 


P 


1 


t 


p 


CP 


rv* 


(V 


• 


» 


p 


p 


i 


a 

© 

CD P 


CP CP 

w 

OJ CVJ 


OJ Pd 

• 

• • 


• *> 

p 

P- P 


P P 


1 1 


1 I 

© 

CO P 


CP CP 

CVJ 

OJ OJ 


CVJ OJ 

* 

• • 


• » 

p 

P P 


p p 


t a 


a t 

© 

© p 


CP CP 

CVJ 

Oi CV 


CV OJ 

• 

• * 


■9 9 

p 

p p 


p p 


ft a 


I 1 

© 

O P 


*-« CP 

CVJ 

cv cv 


r> Cv. 

• 

• • 


* 9 

p 

p p 


P P 


» ft 


1 1 

*-4 

© p 


in <p 

OJ 

P CV 


P (V! 

• 

• • 


9 9 

p 

<> P 


© P 


i 


1 

H 

p p 


in cp 

CVI 

4 w 


■4- OJ 

• 

• 9 


9 9 

P 

< p 


* P 


m a 


cn a 


ro p 


^<p 

rv> 

© cv 


© CVJ 

• 

• • 


9 9 

p 

n p 


fO P 


in a 


in a 

rH 

m p 


0* 

fu 

CD CVJ 


CO OJ 

• 

• • 


• 9 

p 

»^P 


»“• p 


in a 


in a 


V) 


in 


in 

>- 

cn 


IU 

p 

UJ 


CL 


<x 


P 

2 

p 


« 

M 

in 



Ll 




2 




M 




u 




p 



99 



EXAMPLES OF GRAPHICAL COMPUTER OUTPUT PLOTS 


This section of the appendix presents computer-plotted figures of 
typical results calculated by the program for six sets of welding parameters 
shown in Table III (Section IV) . Each figure has pertinent information listed 
to describe the conditions used for a particular calculation. All calculations 
are for 2014-T6 aluminum plates of infinite length and an effective plate 
width of 8 inches or the half width of 4 inches. 

The first 14 figures use the same conditions as specified by NASA for 
their experiments, and thus theoretically simulate these same conditions. 

The first 10 figures (5 pairs) illustrate the temperature and stress 
distributions across the plate at the instantaneous time given on each pair of 
figures. Figures A- 1, A-3, A-5, A-7, and A-9 show the temperature distri- 
butions at 7, 9, 15, 50, and 1000 seconds, respectively. Figures A-2, A-4, 
A-6, A-8, and A- 10 show the stress distributions at these instances. For this 
sequence of figures the heat source is located at the 9-second time position. 

Figures A- 11, A- 15, A- 19, A-23, A-27, and A- 31 illustrate the temper- 
ature (of above ambient) distributions for six sets of conditions in the x- 
direction along planes down the center line of the weld, 1 inch out from and 
parallel to the center line, and down the edge of the plate parallel to the 
center line. Figures A- 12, A- 16, A-20, A-24, A-28, and A-32 illustrate the 
corresponding stress (ksi) distributions'. 

Figures A- 13, A- 17, A-21, A-2 5., A-29, and A- 3 3 plot temperature 
(°F) distributions as they appear out across the plate. 

Figures A- 14, A- 18, A- 22, A-26, A-30, and A- 34 plot corresponding 
stress (ksi) distributions. These contours are essentially the isotherms and 
isostress lines of the plate assuming constant temperature through the thick- 
ness of the plate. These stress contours are unique to this study and provide 
an enlightening picture. 
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RLUniNlM 2014 T6 


FIGURE A- 1. TEMPERATURE DISTRIBUTION ALONG THE SECTION 
AT T = 7, CONDITION 1 
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FIGURE A-2. STRESS DISTRIBUTION ALONG THE SECTION AT 
T = 7, CONDITION 1 
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FIGURE A-7. TEMPERATURE DISTRIBUTION ALONG THE SECTION AT 
T= 50, CONDITION 1 
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FIGURE A-9. TEMPERATURE DISTRIBUTION ALONG THE SECTION AT 
T = 1000, CONDITION 1 
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FIGURE A- 10. STRESS DISTRIBUTION ALONG THE SECTION AT 
T = 1000, CONDITION 1 
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50 SECONDS 



FIGURE A- 12. STRESS CHANGES ALONG THREE LONGITUDINAL LINES 
CONDITION 1 
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FIGURE A- 13. ISOTHERM PATTERN, CONDITION 1 
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FIGURE A- 14. ISOSTRESS PATTERN, CONDITION 1 
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FIGURE A- 15. TEMPERATURE CHANGES ALONG THREE LONGITUDINAL LINES 
CONDITION 2 
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FIGURE A- 18. ISOSTRESS PATTERN, CONDITION 2 
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FIGURE A- 19. TEMPERATURE CHANGES ALONG THREE LONGITUDINAL LIN ES 
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FIGURE A-21. ISOTHERM PATTERN, CONDITIONS 
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FIGURE A-22. ISOSTRESS PATTERN, CONDITION 3 
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FIGURE A-23. TEMPERATURE CHANGES ALONG THREE LONGITUDINAL LINES 
CONDITION 4 
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FIGURE A- 24. STRESS CHANGES ALONG THREE LONGITUDINAL LINES 
CONDITION 4 
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FIGURE A-26. ISOSTRESS PATTERN, CONDITION 4 



FIGURE A-27. TEMPERATURE CHANGES ALONG THREE LONGITUDINAL LINES 
CONDITION 5 
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FIGURE A-28. STRESS CHANGES ALONG THREE LONGITUDINAL LINES 


SO SECONDS 
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FIGURE A- 30. ISOSTRESS PATTERN, CONDITION 5 
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FIGURE A-31. TEMPERATURE CHANGES ALONG THREE LONGITUDINAL LINES 
CONDITION 6 
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FIGURE A- 32. STRESS CHANGES ALONG' THREE LONGITUDINAL LINES 
CONDITION 6 
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FIGURE A-34. ISOSTRESS PATTERN, CONDITION 6 
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